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ABSTRACT 

^^easurements were made o£ tlie limiting shear stress for ttvo naphthenic 
oils of differing molecular weight and tliree blends of the lower 
molecular weight oil and polyalkylmetliaciylate pol)Tners of differing 
molecular weight. The two base oils readied the same limiting shear 
stress for die same temperature and pressure. ITiis was also true for 
all the polymer solutions although the polymer reduced the limiting 
shear stress by about fifteen percent. It seems that limiting stress 
is more a function of material type than viscosity or molecular weight. 

A new falling body viscometer was constructed to operate to 230C and 
to 0.6 GPa, Anotlier viscometer was constructed to extend tlie pressure 
range to 1.1 GPa. 

A new concentrated contact simulator has been developed which allows 
recording of the traction force while the slide-roll ratio is con- 
tinuously varied and the rolling speed is maintained essentially 
constant by a single drive motor. The configuration is that of a 
croivned roller against a disk. 

Measurement of lubricant minimum film thickness of elliptical EHD 
contacts of various aspect ratios were made by optical interferometry. 

ITie data collected wei’e used to evaluate the Hamrock and Dowson minimum 
film thickness model over a range of contact ellipticity ratio where 
the major axis of the contact ellipse was aligned both parallel and 
perpendicular to the direction of motion. A statistical analysis of 
the measured film thickness data showed that on the average tiie experi- 
mental data were thirty percent greater than the film thickness predicted 
by the model. 

Preliminary development of the application of a scanning infrared radi- 
ation system to a tribo-system was completed. A commercial scanning IR 
detector was anployed to measure tlie radiation from an especially con- 
structed elastohydrodynamic device and from a commercial Timken type 
tester. The scanning rate was twenty -five frames per second. 

An analytical study was undertaken for the film thickness developed 
in line contacts under pure rolling conditions employing the limiting 
shear stress rheological model in the various lubrication regimes 
(rigid surface-isoviscous lubricant, rigid surface -variable viscosity 
lubricant, elastically deforming surface-isoviscous lubricant, and 
elastically deforming surface -variable viscosity lubricant). Compari- 
sons were made in all cases with the corresponding Newtonian lubricant 
cases. 
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I. INTRODUCTION 

In previous reports [1] it has been shown that at high pressure, typical 
lubricants undergo a viscoelastic transition from liquid to solidlike 
behavior. Tlie transition pressure was determined for various temperatures 
and was shown to be a function of tlie characteristic rate of change of 
the environment. The constant I’ate transition was shown to occur at 
essentially constant viscosity by plotting isoviscous curves on a pressure- 
temperature transition diagram. 

Investigations f 2 ] of the ultimate shear strength of solidified lubricants 
and of tlie viscosity of the liquids at pressures just below their transi- 
tion formed the basis of a rheological model employing three primary 
lubricant properties; 1cm shear stress viscosity, limiting elastic shear 
modulus, and the limiting shear stress. The first of these properties 
is familiar and has been reported often in the literature. The limiting 
elastic shear modulus is less familiar but has been measured for many 
years by several techniques including ultrasonics. Tiie limiting shear 
stress of liquid lubricants has been the subject of speculation for many 
years and indeed the nature of traction in EID) contacts has led several 
researdiers to support that view for some time. Conparison of the model 
with traction versus slide-roll ratio data from other researchers [ 3 ] 
and with traction versus sideslip data from our laboratory showed that 
the model was essentially correct for isothermal conditions. 

Since tlie lubricant was usually in the aiprphoas solid condition when 
the limiting shear stress was measured, it was e:^ected that the history 
of the formation of tlie solid influenced the limiting stress values 
measured. In tl.e previous report fl] it was shoivn that the isothermal 
compression history produced a limiting stress lower (and more charac- 
teristic of EIID) tlian the previously used isobaric cooling process. The 
pressure range was extended for limiting stress and shear modul^ measure- 
ments and the capability of measuring these properties with solid polymers 
was added. The yield shear stress behavior of the solid polymers was 
found to be conparable to that of the liquid lubricants. 

Daring the current year measurements were made of the limiting shear 
stress for two naphthenic oils of differing molecular weight and three 
blends of the lower molecular weight oil and polyalkylmethacrylate 
polymers of differing molecular weight. Tlie two base oils reached the 
same limiting shear stress for the same temperature and pressure. This 
was also true for all the polymer solutions although the polymer reduced 
the limiting shear stress by about fifteen percent. It seems that 
limiting stress is more a function of material type than viscosity or 
molecular weight. 

'I A new falling body viscometer was constructed to operate to 230C and to 
0.6 GPa. Another viscometer was constructed to extend the pressure range 
to 1.1 GPa. 
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• A new concentrated contact simulator has been developed which allows 
recording of the traction force while the slide- roll ratio is con- 
tinuously varied and the rolling speed is maintained essentially con- 
stant by a single drive motor. The configuration is that of a croimed 
roller against a disk. Preliminary results are promising and further 
development is expected to yield the capability of measuring the effect 
of kinematics and lubricant properties on the linear, low slide roll 
ratio portion of the traction curve. 

' Measurement of lutiricant minimum film thickness of elliptical EHD 
contacts of various aspect ratios were- made by optical interferometry. 

The data collected were used to evaluate the Hanrock and Dowson minimum 
film tliickness moclel over a range of contact ellipticity ratio where 
the major axis of the contact ellipse was aligned both parallel and 
perpendicular to the direction of motion. A statistical analysis of 
the measured film thickness data shoived that on the average the experi- 
mental data were thirty percent greater than the film thickness predicted 
by the model. 

' Preliminary development of the application of a scanning infrared radia- 
tion system to a tvibo-system was completed,/ A commercial scanning TR 
detector was enployed to measure the radiation from an especially con- 
structed elastohydrodynamic device and from a comnercial Timken t)Tpe 
tester. Tlie scanning rate was twenty-five frames per second. The effort 
described in this report obtained data by photographing the CRT display. 
Currently, efforts are underway to record the video signal for post analysis 
which is expected to greatly enliance both the spatial and time resolution 
of the information. 

An analytical study was undertaken for the film thickness developed in 
line contacts under pure rolling conditions employing the limiting shear 
stress rheological model in the various lubrication regimes (rigid surface- 
isoviscous lubricant, rigid surface- variable viscosity lubricant, elastically 
deforming surface-isoviscous lubricant, and elastically defoiming surface- 
variable viscosity lubricant) . Conparisons were made in all cases with 
tlie corresponding Newtonian lubricant cases. 
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II . LUBRICANT SHEAR RHEOLOGICAL MEASUREMENTS 
Experimental Apparatus 

Lubricant shear rheological measurements were conducted on several 
materials and in several devices. The devices were botli those for 
measuring primary properties (viscosity and limiting shear stress) 
and a traction rig for measuring traction versus contact kinematics. 

Some of the devices were constructed as part of this year’s effort. 

Several pieces of apparatus were employed. They consist of devices 
to measure limiting shear stress, viscosity and traction. The latter 
two, as explained, were constructed as part of this year’s effort. 

Measurements of shear properties of fluids in tlie amorphous solid 
condition were conducted in the Constant -Pressure Stress -Strain 
Apparatus which is described in the previous year’s annual report [1]. 

A new falling body viscometer (Figure 1) was constructed to operate 
to 230C and 0.6 (Pa. The higher tenperature capability is due to new 
designs of the linear variable differential transformer and the vessel 
closure seal. In addition, the falling body now translates within a 
loose fitting sleeve inside the pressure vessel so that changes in 
vessel bore diameter with taiperature and pressure do not affect the 
viscometer calibration. The sleeve is closed at one end and sealed 
at the other by an isolating piston forming a removable viscometer 
cartridge which holds about 2 cm^ of sample. 

In addition, a new viscometer was constructed to extend the pressure 
range to 1.1 GPa. In this design (shown in Figure 2) the pressure 
generating intensifier is incorporated in the viscometer and supports 
a portion of the outside diameter of the soft non-magnetic vessel for 
use above lOOC the intensifier is forced air cooled. A viscometer 
cartridge similar to that described above is used. 

A new concentrated contact simulator (Figure 3) was developed and is 
nearing completion which allows recording of traction force and infrared 
surface temperature while the slide-roll ratio is continuously varied 
and the rolling speed is maintained essentially constant by a single 
variable speed drive motor. The configuration is that of a crowned 
roller against the face of a disc. The disc is loaded against the 
roller by weights acting through a thrust bearing and keyed to prevent 
rotation. The sideslip angle is adjusted from approximately zero to 
0.05 radians by means of a micrometer which tilts the disc rotational 
axis about a lower sipport bearing. The angular velocity of the roller 
and disc are coupled at a fixed ratio while the sliding velocity is 
varied by varying the distance from the contact to the axis of the 
disk. A change in slide- to-roll ratio (Z) of 0.14 can be continuously 
scanned about a mean value of slide- to- roll ratio which is variable in 
larger increments. A preliminary traction curve is presented in 
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Figure 4 wluch shows the effect of side slip angle. 


Rl\eo logical Exi^eriinents and Results 

Flow curves (shear stress vs. shear rate) were generated using tlie 
constant pressure stress -strain apparatus for two naphthenic mineral 
oils of differing molecular weight and viscosity and three blends of 
tlie Im viscosity oil and four percent by weight of polyalkylipthac- 
lylate pol>Tners. Hie polymer -molecular weights were 3.5 x 10^ , 

6 X 10^, and 20 x 10^. The material history for measurements in tlie 
amorphous solid regime was that of nearly isothermal compression 
from the liquid phase. Fluid descriptions can be found in Appendix A. 

In Figure 5, two flow cliarts are presented for tlie low viscosity 
naphtlienic base oil (R620-15) for pressures of 1.09 and 0.93 GPa 
respectively. In eadi case the temperature is varied. At the higher 
temperatures linear viscous behavior is evident at the lowest shear 
rates. Also, increasing the tenperature apparently lowers the limit- 
ing shear stress. Figure 6 shows the shear behavior of both the low 
viscosity base oil (R620-15) and the high viscosity base oil (R620-16^ 
at 26C and various pressures. Included are results for tlie low vis- 
cosity oil at 46C. It can be seen that lowering tlie pressure reduces 
the lijidting shear stress and brings on viscous behavior at lav shear 
rate. Fluid descriptions can be found in Appendix A. 

Blends of polyalkylmethacrylate and the low viscosity oil were prepared 
and run in the stress-strain apparatus at 26C (Figure 7). Aside from 
viscosity differences tlieir behavior is apparently tlie same. 

In Figure 8 curves are drawn exclusive of the data points for the two 
base oils and three polymer blends at 26C and .93 (^a. For this par- 
ticular temperature and pressure, the base oils readied essentially 
the same limiting stress. This was also true of the polymer solutions 
although the polymer reduced tlie limitin'' stress of tlie base oil by 
about fifteen percent. It seems tliat the limiting shear stress is 
more determined by material type than by viscosity or molecular weight. 

Following conpletion of tlie high temperature viscometer, a pressure 
viscosity isotlierm was generated for 5P4E (polyphen)'’! ether) at 227C 
(440F) to 0.5 GPa as shown in Figure 9 along with previously presented 
low temperature data. In Figure 10 are four isotherms for R620-15 
developed witli the instruments indicated on the figure. At 149C data 
was taken to a pressure of 1.07 GPa. 

Figure 11 is a collection of pressure-viscosity isotherms for the two 
base oils and two of the polymer blends at 29, 99 and 227C. Pressure 
viscosity coefficients for these materials are tabulated in Table 1. 
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Since tJie lav shear rate portion of the flav cliarts in Figures 4, 5 
ami 6 shows viscous heliavior, it is informative to plot tlie low stress 
viscosit)' calculated fran these cliarts on tlie respective isothems 
shown in Figure 10. This was done in Figure 12. Altliough there is a 
large pressure raiige in whicli no data is available, the extrapolated 
isotliems from the falling bod>' exiiorimonts can be reasonably extraiio- 
lated to fit the data from tlie stress-strain apparatus. This provides 
us with a fin’ther tool foi' extending tlic limits on pressure and vis- 
cosity with respect to viscosit)' measurement. Note that at veiy high 
pressure tlie poUnner thickener has reduced the viscosity of the naph- 
thenic solvent, although at lav pressuix? it has increased it. This 
unexiiected result is seen to have been oliserved with botli instruments. 
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III. AN EXPERIMENTAL EVALUATION OF THE HANCOCK 
AND DaVSON MINIMUM FILM THICKNESS EQUATION 
FOR FULLY FLOODED EHD POINT CONTACTS 


Introducti on 

Tlie researdi reported in lliis sectuM. was the result of tlie Masters 
degree tliesis researdi of Mr, K, A. 1‘oye [4] a.id vdll be presented at 
the ASME/ASLE Intemational Lubric.it ton Confeicnce in August 1980 
and published in Transactions ASME, Journal of Lubrication Tedinol- 
ogy in substantially tlie form presented here. 

The minimum thickness of the lubricating film separating the surfaces 
in an EHD contact is of primary interest to a design engineer con- 
cerned with sudi contacts. The successful separation of the surfaces 
under the anticipated operating conditions can significantly increase 
tlie useful life of the madiine parts involved. If separated, tlie 
failure mode will be tlie long term limitation of surface fatigue. 

If the surfaces are not separated adequately, surface interaction 
will shorten the contact life and may even result in sudden scuffing 
failure. 


The puriiose of tliis study was to determine experimentally, for a prac- 
tical range of elliptic! ty ratio, the effect of contact geometry on 
minimum film thickness in a fully flooded isotliermal EHD point contact 
in pure rolling. Further, it is intended to use tlie data collected to 
statistically examine, by regression analysis, tlie Hamrock and Davson 
minimum film tliickness equation 



3.63 ij0.68j.0.49 ^^-0.073^ _ g-0.68k^ 


( 1 ) 


Hamrock and Dows on [ .S-8 ] only intend tliis equation be used for 
ellipticity ratios greater than one (k ^ 1) and it is speculation on 
the authors' part tliat it is valid for ellipticity ratios less than 
one (k < 1) . As tlie experimental results sliow , this speculation is 
justified. 


In addition, tlie applicability of the Hamrock and Dowson film thick- 
ness equation to contacts aligned witli tlie major axis in the direc- 
tion of rolling is evaluated. 

Tlie objective of the experimental work performed was to obtain EHD 
film tliickness measurements for a number of ellipticity ratios , for 
eadi of tlie above alignments of tlie major axis, at various levels of 
dimensionless speed and load. The dimensionless material parameter 
was maintained constant. Tlie dependent variable, dimensionless film 
thickness, was obtained for eacli set of independent variables: speed, 
load, and ellipticity ratio (or radius ratio), £ind tlie results then 
compared witli the predictions of the Hamrock and Dowson model. 
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Nomenclature 


a contact ellipse dimension perpendicular 

to the direction of motion 


k = a/b 




l-v2 l-v2 

“ET* Trr 


F 

G 5 a E’ 

y u 

u = _2 

^ " E’R . 






*Wn.exp. 


min 


contact ellipse dimension parallel to 
the direction of motion 

ellipticity ratio 

Ej^, E^ and are elastic moduli 

and Poisson’s ratio for body 1 and 2 

load 

dimensionless material parameter 
dimensionless speed pai'ameter 

dimensionless load parameter 

film thickness 

measured minimum film thickness 
minimum film thickness calculated by 



Iv 


min 



dimensionless minimum film thickness 


r ,r roller radii parallel (x) and 

^ perpendicular (y) to direction of motion 

R ,R equivalent principle radii of curvature 

^ of the contacting surfaces parallel (x) 

and perpendicular (y) to the direction 
of motion 
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■v* 



where 


^2' V’ V 

u s I (Uj + U 2 ) 


are the radii of surfaces 1 and 2 in 
the direction of motion (x) and per- 
pendicular to it (y) 

average rolling velocity 


u 


1»“2 


velocity of surface 1 and 2 


absolute viscosity at HID inlot temperature 


a 


pressure viscosity coefficient at inlet 
temperature 


j 
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Experiniental Apparatu s 

EHD Simulator 

The apparat 3 used to create conditions similar to those encountered 
by tlie contacting surfaces in rolling element bearings and gear teeth 
is shown in Figure 13. Tlie principal parts of tlie EHD simulator are 
the contacting surfaces of a hardened steel crojvned roller and a flat 
syntlietic sapphire disk. 

Hie croitfned roller is mounted in pillow blocks containing radial ball 
bearings (9.525 mm I.D.) on an aluminum block with a slot milled just 
below the surface of tlie cravned roller whidi serves as an oil batli. 
The crowned roller support is shown in Figure 14. 

The sapphire disk is mounted in a ball thrust bearing (11.0 cm O.D.) 
and is driven about an axis perpendicular to that of tlie roller rota- 
tion. Hie sapphire is driven through a bevel gear drive by a d.c. 
servo motor (0 . 15 HP) . 

Hie sapphire in its thrust bearing and its complete drive system are 
mounted on tlie contact loading lever. Hie load at tlie contact was 
4.4 times the load at the end of the lever arm. A viscous dancing 
mechanism was used to eliminate meclianical vibrations and keep the 
load constant. 

The speed of the sapplare was measured by an optical taclicaneter which 
was connected to the drive motor shaft through a flexible coupling. 

The tachometer provided an input signal of fifty counts per revolution 
to a digital counter. 

To insure near pure rolling conditions (surface velocities equal) tlie 
speed of the crowned roller was also monitored. Pure rolling was 
desired to avoid the increase in lubricant temperature and subsequent 
decrease in inlet viscosity due to the viscous energy dissipation 
associated with sliding. A bifurcated fiber optics bundle was aimed 
at an extension of the roller shaft to ivhich two bars of reflective 
plastic were glued 180° apart. The light transmitted and received 
by the fiber optics bundle was conditioned into an electronic signal 
of two counts per revolution as input to a digital counter. Hie bulk 
fluid tenperature was also monitored witli a tliennocouple attaclied to 
a digital tliermometer. 


Description of the Bearing Surfaces 

Hie crowned rollers were cut from 3.81 cm diameter A-2 tool steel bar 
stock. A-2 tool steel was selected for its machinability and harden- 
ability. In seme cases cutting tools with the desired crown radius 
were milled from 0-1 tool steel and then used to cut the radius into 
the roller surface. For tlie smaller crown radii, less than 7.9 mm, 
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and the larger, greater than 15,9 mm, a step was cut into tlie roller 
of approximately the desired crown diameter in widtli. This step was 
then sanded whicli resulted in a curved surface. The finished rollers 
were then hardened to a Rockwell hardness of Rc 60. The hardened 
rollers wei’e then sanded and polished with emery and polishing paper 
to a surface roughness of less than 0.13 \m (Arithmetic Average). 

The actual finished crown radius of a roller was determined by mea- 
suring the ellipticity ratio of the elliptical contact produced by 
loading the sapphire disk against it in the EUD simulator. The rolling 
radius of the roller was measured with a micrometer. Thus knowing the 
ellipticity ratio and the radius in the direction of motion the crown 
radius could be determined from the relationship , 

k = 1.0 339 (R^/R2)°‘^^^, ^1-^2 ^ ^ 

derived nun.erically by Brave and Hamrock [ 9 ] from Hertzian contact 
stress theory. 

The geometry of the crowned rollers used in tliis investigation is given 
in Table 2. The physical properties of A- 2 tool steel are given in 
Appendix B. 

The synthetic sapphire (Al20j) disk was chosen for its strength and 
hardness. In addition, it has the optical properties, transparency 
and refractive index, necessary for optical film thidoiess measure- 
ment. The disk measured 8.89 an diameter by 0.32 cm thick with a 
surface rouglmess of 0.00635 vim A. A. and optically flat to within an 
eighth of a wavelength. The EHD contact surface of tlie disk was 
coated with Inconel to give a partially reflective surface. 

The lubricant used was a naphthenic base mineral oil (Sunoco R-620-16). 
It was selected for its base viscosity at room temperature, whidi 
would insure an adequate film formation in the smaller contacts. 

The properties of the lubricant crucial to the film thickness equa- 
tion, namely, inlet viscosity and pres sure -viscosity coefficient, 
were measured in the Tribology Laboratory at Georgia Tech at tem- 
peratures near the operating temperature. The results of these mea- 
surements are given in Appendix B, along with otlier physical proper- 
ties of the lubricant. 


Qjtical Interferonetry Film Thickness Measurement Teclmique 

The optical interferometry technique for measuring EHD film thickness 
which was used in tliis investigation was developed by Foord, Wedeven, 
Westlake, and Cameron [10]. Tlie technique was adapted for use on the 
EHD simulator used by Nagaraj [11]. 
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A didiromatic liglit source was directed onto tlie vertical illuminator 
of a Leitz Metallurgical microscope to whidi a Polaroid camera had 
been mounted. The film thickness in the contact could be deduced 
from a color photograph of the interference pattern of tlte dichro- 
matic light reflected up into tlie camera from the contact. 

The optical arrangement used in the experimental investigation is 
very similar to that of Foord et al. [10 ] . A dichromatic light source 
was used whicli consisted of a 5460 A green band and a 6300 A red band. 
Tlie dichromatic incident light source produces an interference pattern 
which consists of three predominant bright fringe colors and one dark. 
The three bright fringes were red, green and yellow, and the dark 
fringe appeared black. In between the sharp, distinct, bright and 
dark fringes there appeared shades of the predominant colors. The 
sequence in whicli these colors appeared varied over 2.3 pm of film 
thickness and then repeated. 

The effect of tenperature and pressure on tlie refractive index of the 
lubricant has been considered by Nagaraj [11]. Nagaraj reported tlie 
variation in refractive index due to pressure variations over the 
contact area to be less than ten percent for a system very similar 
to tlie one used in this investigation. In addition he reported a 
two percent change in film thickness for a 50 percent error in Hertz 
pressure estimate at 0.65 GPa, and a 1.5 percent change in refractive 
index for a temperature rise from 38C to lOOC at 0.138 GPa. There- 
fore, the effect of temperature and pressure on refractive index of 
the lubricant were neglected in the calculation of film thickness. 

The optical system used in the experimental measurements of EHD film 
thickness is slioivn schematically in Figure 15. A continuous collimated 
light source directs a beam of light down an optical bench into the 
vertical illuminator of tlie microscope. An aperture was used to 
control the amount of light from the 30 watt, variable intensity light 
source. The light then passes through a dichromatic filter and a 
condensing lens system. 

The continuous light source was used to focus the interference 
pattern of tlie static EHD contact and observe films of contacts 
where the surface velocity was less than 10 cm/s. The EHD contact 
became increasingly more dynamic in the field of view of tlie micro- 
scope witli increasing rolling speeds. This motion of the EHD contact 
was due to slight mechanical misaligranents which are magnified by tlie 
microscope. At speeds over 10 an/s it was necessary to photograph 
the interference pattern reflected from the contact with an exposure 
time short enough to stop the motion of the contact. A flash unit 
with flash durations from 30 to 40 microseconds was found to provide 
the correct exposure for Polaroid Type 58 Land Film when the flash 
unit’s thyristor circuitry control for flash duration was set at 
manual. A Polaroid four indi by five incli film holder and camera 
were mounted on top of the microscope. 
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Experimental Teclmique and Method of Analysis 
Design of the Exi^eriment 

The experimental procedure used to obtain minimum film thickness 
measurements was designed around two primal^' considerations, These 
considerations were first to use a practical range of dimensional 
speed, dimensional load, and ellipticity ratio in obtaining film 
thickness data from Effl) contacts, and second, to use a sufficient 
number of data points within eadi range to obtain a meaiiingful 
regression analysis. 

Tlie Hajurock and Do\\'Son equation for minimum film thickness , Equation 
(13 contains four independent dimensionless variables. Tlie equation 
was developed for values of ellipticity parameter (k) greater tlian 
one. The dimensionless material parameter G was not varied in this 
investigation and was 10451. Thiis, for eacli ellipticity ratio 
examined only tlie dimensionless speed and dimensionless load were 
varied and this variation was achieved primarily by Varying the 
applied load and the rolling velocit)'^ of the bearing surfaces. 

(There was some variation in the rolling radius % caused by finish- 
ing tlie croivned contact surfaces of tlie steel rollers.) 

The Hamrock and Dowson minimum film thickness equation predicts the 
effect of contact geometry on minimum film thidmess in HID contacts 
to be most prominent for tlie range of ellipticity ratio between 0 
and 4. The relationship between ellipticity ratio and minimum film 
thickness as expressed by the Hamrock and Dowson equation is illus- 
trated in Figure 16, which shows less than a ten percent effect of 
ellipticity ratio on minimum film thickness above k = 4. 

Early experimental work concerned witli film thickness in EHD point 
contacts by Arcliard and Cowking [12] reports similar conclusions 
with respect to the effect of contact geometry (particularly side 
leakage) on film thickness in EHD contacts with ellipticity ratios 
in tlie range 0.48 < k < 5.02. Recent work by Baliadoran and Gohar 
[13] shows no significant effect of ellipticity ratio on film thick- 
ness in the range 3.5 < k < 7.5, wliidi sipports the ipper part of 
the curve in Figure 16. Gledliill, Jackson, and Cameron [14] inves- 
tigated low ellipticity ratios (k r 0.1) while studying tlie EHD of 
asperities. Therefore, the range of interest for ellipticity ratio 
in this investigation was chosen as being 1/8 < k < 4. Beginning 
at the low end of this range and moving up, the practical applica- 
tions of this range of ellipticity ratio include: micro-elastohydro- 

dynamic lubrication of surface asperities , worm gear contacts , 
roller-bearing rolling element- race contacts, and contacts in gear 
teeth. 

Dimensional speed and load ranges were diosen to cause the dimen- 
sionless speeds and loads to lie within the ranges used by Hamrock 
and Dowson in the numerical procedure used to develop tlieir minimum 


film thickness equation. These ranges of dir^nsionless speed and 
load were maintained for all but a few cases at the extreme ends 
of the elliptic! ty ratio range. The few cases which did lie out- 
side the desired range of dimensionless speed and load used 
Hamrock and Dowson never did so by more than an order of magnitude. 

The range of dimensionless speed maintained was 

2.14 X 10"^^ < U - ^ < 8.90 X 10"^^ , 

which, for the mechanical system used in this investigation (with 
an average rolling radius taken as R,, « 1.5875 cm) correspond to a 
speed range of, 31 cm/s < u < 128 an/s. The range of dimensionless 
load maintained was 

0.038 X 10'^ < W - S-j < 5.32 X 10'® , 

E-^ 


which corresponded to a range on dimensional load of 2.74 N < F < 

384 N. 

For each ellipticity ratio examined, three levels of speed and three 
levels of load were selected from within the ranges mentioned above. 
This resulted in nine cases of speed anC load for each ellipticity 
ratio, except for a few sets of speed and load which were prohibited 
by equipment limitations . The levels at which speed and load were 
set were selected in equal logarithmic intervals of each other. 

The rolling velocities used were 50, 90 and 125 cm/s for every 
ellipticity ratio at every load except for k « 0.117, where 30, 50 
and 90 cm/s was used. (Photographing this contact was difficult at 
high velocities.) These rolling velocities were selected as being 
at a level of practical applicat.ion arid were well suited to the 
mechanical drive system used to drive the sapphires rotations. 

The three levels of load used for each ellipticity ratio were chosen 
to yield a maximum Hertz pressure in the EHD point contact of 0.69, 
1.03 and 1.55 GPa (100 x 10^, 150 x 10^, and 225 x 10^ psi). 

The loads and speeds used were checked to insure the operating regime 
was in the elastic- variable viscosity (EHD) regime according to the 
elliptical contact regime theory of Hamrock and Dowson [15 ] for 
those cases of k = 1.0, 3.0 for which regime charts are available. 

For the cases of k < 1, it was assumed that Hertz pressures of 0.69 
GPa (100 kpsi) or more should insure operation in the EHD regime. 


Experimental Procedure 

Individual film thickness measurements were made by a set procedure. 
After the load to the contact was zeroed, tiie desired load to the 
contact was applied by placing the appropriate dead weights on the 
loading platform hanging from the end of the loading lever beam. 

The static interference pattern was located and focussed in the 
field of vie\sr of the microscope using the continuous light source. 

The flash unit was then mounted in its position on the optical bench. 

Motion of the contacting surfaces was initiated by sipplying power 
to the d.c. motor. The power was increased until the desired speed 
was obtained as evidenced by the readout on the display of tlie digital 
counter. The rotational speed of the sapphire was checked against 
tliat of the crowned roller to insure they were equal, which indicates 
pure rolling. 

Wlien the desired rotational speeds of the contacting surfaces were 
achieved the resulting interference pattern was photographed by open- 
ing the camera’s shutter and triggering flash unit. The entire 
process , from the moment motion was initiated to the moment the 
photograph was taken, took 60 seconds on ti'e average. The minimum 
film thickness was then determined frtmi the photograph. 

In all, seven crowned rollers with seven different radius ratios 
were used at three levels of speed for each of three levels of load. 
Sixty-three film thickness measurements were individually atten^ted 
by the cdjove procedure for the sixty-three combinations of the three 
independent variables of interest. Fifty-seven useful film thickness 
measurements were obtained. The six cases which did not prove fruit- 
ful were because of extremes in tlie load necessary to generate desired 
Hertz pressure. 


Statistical Analysis of the Experimental Data 

After the minimum film tiiickness for each of the fifty-seven combina- 
tions of speed, load, and ellipticity ratio were obtained, the film 
thickness results and the corresponding independent variables from 
Equation (1) were entered on a conputer data file. This data file 
was used as an iiput file. 

The minimum film thickness was calculated for an EHD point contact 
using the Hamrock and D.xvson equation (Equation 1) given the necessary 
physical property data .and the operating parameters (speed and load) 
for the solid- lubricant-solid system. Also calculated was the 
absolute value of the difference betiveen the measured minimum film 
thickness, h^ , and the calculated value, expressed as 
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and the ratio expressed as 

^mln.exp^inin.^ ’ 

The variables of Equation (1) were transformed into the natural logaritlm 
of these variables for use in a log-multilinear regression analysis. 

The data stored on the output file were used as input for a Biomedical 
Canputer Programs [16] P-series regression routine called PIR. M)P1R 
estimates a multiple linear regression tliat relates a dependent variable 
to several independent variables. Tlie routine uses a least squares metliod 
to estimate the coefficients Sj , 0^ , • • • » 0„ and the intercept a in tlie 
equation x 4 P 



Y « ot + 

* "■ * Vp * ^ 

C2) 

vdiere 

Y - 

dependent variable 



X - 

independent variable p 



e - 

error 



The estimates of the intercept a and the coefficients 0, denoted a and 
b, are found by minimizing the expression 


J (Y - a - bpXp)2 (3) 

for all cases of dependent and independent variables supplied. 

A multilinear regression analysis was performed by BMDPIR with tlie data 
supplied from the output file for the constants a and 0p in the equation 

+ Bj «n(l - e"*’-®®'') . (4) 

Perfect correlation of the experimental data with the Hamrock and 
Dowson model v/ould result in values of 

a » An(3.63) + 0.49 tn(10454) = 5.82 


0.68 
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$2 • " 0.073 

“ 1.0 


BMDPIR qlso performs univariate statistics on all tlie variables 
siq?plied. The statistics calculated include: mean, standard 
deviation, and variance. Tlie coefficients of correlation for tlie 
regression are also supplied. 


Results of the Experimental Investigation 
and Statistical Analysis of the Data 

Experimental Results 

The measured dimensionless minimum film thickness results and the 
corresponding dimensionless speeds and loads for the seven ellip- 
tic! ty ratios examined are presented in Tables 3 through 9 . In 
addition, the absolute value of the difference between the measured 
dimensionless minimum film thickness and the value predicted by the 
Hamrock and Dows on Equation (Equation 1) , along with the quotient 
of the two, are also presented. Exanples of the photo -micrographs 
used to measure minimum film thickness for each of the seven ellip- 
ticity ratios examined are presented in Figure 17 through 23 corre- 
sponding to rollers in Tables 3 through 9 respectively. The 
original photomicrographs were in color. The reproductions shown 
are from black and white high contrast photographs of the original 
with unfiltered light. In Figures 17 and 18 the flow is from bottom 
to top while in the remaining figures it As from right to left. In 
all cases the minimum film thickness is on the side or trailing edge 
depending on conditions. 

Figure 24 is a plot of measured dimensionless minimum film thickness 
versus the dimensionless film thickness predicted by the Hamrock and 
Dowson model. Tlie line H . = represents perfect agreement 

between the experimental results anfl®Tne predictions of the Hamrock 
and Dowson model. 


Results of the Statistical Analysis 

The multilinear regression estimated the coefficients a, B-j^, &2 
Bj in the equation 

toCHj^n.exp^ = a + 3i An (U) + 83 An(W3 + BjAnCl - e'°‘^®^)(5) 
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a * 7.344 

3^ 0.751 ± 0.045 

32 = - 0.097 t 0.021 

3j = 0.955 * 0.010 


Thus, the tnodel developed by the multilinear regression analysis was 


^min “ u0.751 q0.49j^,- 0.097 _ g-0.68k^0.955 ^ 

2 

The correlation coefficient, k, product of correlation, R , and the 
standard error of the estimate for the multilinear regression were 


R = 0.9856 
R^ = 0.9713 


STANDARD ERROR OF ESTIMATE = 0.1439 


where a correlation coefficient of one and a standard error of the 
estimate of zero indicate perfect correlation. 

The univariate statistics performed by BMDPIR were of particular 
interest for the variables 

l^n " \iin.exp^ 

and 

^^min. exp^^min^ 

The mean and standard deviation of these two variables are given in 
Table 10. 

The ratio of experimentally measured to calculated film thickness was 
also examined for each ellipticity ratio. The mean and standard devi- 
ation for this statistic is shown in Table 11. As seen from tliese data 
there is somewhat of a trend for the experimental film thicloiess to 
deviate more from the calculated value for the lower ellipticity 
ratios. This is not necessarily attributed to the lack of predictability 
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of the equation at low ellipticity ratio. It can also be attributed 
to ^e fact that tlie measurement teclinique has a minimum uncertainty 
wliidi is a larger fi'action of tlie measured film tliickness as the film 
tliickness decreases, as it does for the lo^^r ellipticity ratio cases. 

The measurement of tlie film tliickness is more difficult at tlie loiver 
ellipticity ratio (light loads, smaller interferograms) also. If we 
restrict the conparison to ellipticit)’^ ratios of appi’oximately one 
or more, as Uamrock and Dov\?son intended for their equation, tlie mean 
ratio of experimental to calculated film tliickness is 1.23 and tlie 
standard deviation about nine percent. Tliis is only slightly better 
than tlie overall statistics of 1.30 and 0.11 resiiectively. 

Also examined in the statistical analysis were the residuals based on 
tlie generated regression equation (6) . These residuals and tlieir 
squares were plotted as functions of the dependent variable (calculated 
film thickness) and eacli of the independent variables (ellipticity 
ratio, sjieed, mid load). There were no obvious trends of the residuals 
with miy of the vjiriables except that the spread of them was larger 
at tlie smallest film thicknesses as might be exiiected. Even in that 
case the residuals wei'e centered aroiuid zero. 


Discussion 

The results presented show a reasonable agreement bet\>rcen tlie measured 
film thickness data and the Mamrock and Dowson model. 

The regression coefficients estimated by the multilinear regression 
analysis did not differ significantly from the corresponding exiionents 
of the Uamrock and Dowson minimum film tliickness equation but tlie 
multiplicative constmit was greater by a factor of 4.5. 

M . , . 3.63u‘’-f’V-«W"-'’”(l - (1, 

luin 

The model developed by the multilinear regression analysis. 


inin 


16.6 u«-751c0.4V0.097^^ 


„-0.68k.0.955 
e ) 


( 6 ) 


fit tlie exjierimental film tliickness data well; the correlation coeffi- 
cient was R = 0.9856. The mean measured minimum film thickness was 
witliin one percent of the multilinear regression model (Equation 6) 
prediction witli a standard deviation of 13.98 percent. Tlie F~ ratio 
from this regression analysis was 666.767 with three degrees of freedom 
in tlie regression and 59 degrees of freedom in tlie residual. Applying 
tlie F- ratio test, this iinjilies that the independent variables used 
(speed, load and ellipticity ratio) account for the variation of the 
dependent variable (film thickness) at tlie 99.999 percent confidence 
lindt (P(tail) = 0.00000). 
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Hvidence of the correlation between the iiieasured film thickness data 
and tlie Hainrock rind Denison model can be found in Table 10 and 11, The 
mean difference between measured and predicted minimum film thickness 
was 6.682 pin. (0.17 phO . Tlie mean percent difference was 30.2 per- 
cent and tlie standard deviation was 21.4 percent. The mean film 
thickness measured was 26.27 pin. (0.667 pm)* procedure for 
measiu'ing the minimum fibn thickness from the photondcrograph of 
the interference patteni is only good witJiin plus or minus 2.5 micro- 
incJies (0.0127 pm). Tlius, for a mean measured film thickness of 26.7 
microindies tlie exiierimental procedure came witliin 30.2 percent of tlie 
Hamrock mid Doivson equations prediction, 8.01 pin, of whicli only 2.5 
microindies could be attributed to lack of accuracy' in measuring the 
film tliickness. The data presented in Table 10 show that on tlie average 
tlie measured film tliickness data was thirty percent more tlian the 
Hamrock and Dowsou predictions and 68 percent (i.e., ± one standard 
deviation) of all of the measured iidnimum film thickness data were 
from ten percent more than to 50 percent more tli.Tn tlie Hamrock and 
Dowson predictions. (Assioning a normal distribution of • 

Figure 24 presents the dimensionless minimum predicted film mii^iess 
from tlie Hamrock and Dowson equation and the corresponding exiierimental 
measurement. It illustrates the fact that the measured miniminu film 
thicknesses wore greater than tliose predicted. 

However, if the conqiarison i.s restricted to the ellipticity range for 
which tlie Hamrock and Dowson equation was intended (k - 1) , it is 
scon from Talile 11 that the agreement betiveen exiieriment and tlieory is 
better ill that the mean of the measured film thickness is only 15 per- 
cent greater than the prediction mid tlie standard derivation is also 
smaller (about eight percent). Therefore in the range of ellipticity 
ratio for whicli the equation was develojied, plus and minus one standard 
derivation of the exiierimental data range from about seven percent to 
about 23 percent more tlian the predicted values. 

As indicated it was speculative to exiiect the Hamrock and Dowson 
equation to be applicable for ellipticity ratios less than one; a 
range which was not considered by Hamrock and Dowson in developing 
tlie equation. However, as seen from liquation (6), the multilinear 
regression model developed from tlie data, tlie exponent on tlie ellip- 
ticity fimction for tlie best fit I'egression is 0.955 when the entire 
ellipticity range is included (0.117 < k f 3.70) compared to an 
exiionent of one on tlieir model. This is particularly significant 
when it is recalled that the ellipticity I'atios in this work were the 
experimentally measured values. 

It is interesting to point ouc that the value of 0.751 for the exponent 
on the dimensionless speed is close to the value of 0.744 arrived at 
exjierimentally by Arcliard for point contacts , and through analytical 
methods by Grubin, Crook and Dowson and Higginson for line contacts [17]. 


20 


The exponent on the load term from the regression model is somewhat 
greater than that in the analytical equation (i.e., 0.097 compared 
to 0,073) which indicates the experimental data exhibited a somewhat 
greater dependence on load than expected. 

In this experimental work every effort was made to assure conditions 
that were assumed in developing the Hamrock and Dowson equation 
except for considering k < 1 and the material parameter (G = 10451) 
which was greater than Hamrock and Dowson range. Of their 34 cases 
considered 31 were for G » 4522 and one each at 2310, 3491, 6785. 
There was a good knowledge of the properties of the lubricant, par- 
ticularly viscosity and pressure viscosity coefficient, at the opera- 
ting conditions. The lubricant was Newtonian. Themal effects were 
minimized to satisfy the isothermal assunption. The roller support 
assembly used was free of any extraneous frictional effects at the 
contact surface. In addition, by photographing the contact inter- 
ference pattern within one minute of the beginning of motion, in- 
dividually for eadi set of conditions examined, little time was 
allowed for bulk thermal effects to develop in a rolling contact, 

(as opposed to a sliding contact) . The bulk fluid tenperature of 
the lubricant was never noted as having varied more than 0.5C (1°F) 
about 23C (74F). In spite of these precautions, the mean measured 
minimum film thickness was 30 percent greater than the predicted 
values calculated from the Hamrock and Dowson equation (1) . 


Conclusions 

The collective results of the evaluation of the Hamrock and Dowson 
model for minimum lubricant film thickness in elastohydrodynamic 
point contacts support the model. In addition, there was no sig- 
nificant deviation of the model from the experimental data for any 
particular range of ellipticity ratio examined. The mean measured 
minimum film thickness was 30 percent greater tlian the Hamrock and 
Dowson prediction (with a standard deviation of 21.4 percent). There- 
fore, it is the conclusion of this evaluation the Hamrock and Dowson 
model conservatively predicts the minimum film thickness in EHD con- 
tacts within the accuracy normally required in mechanical design 
calculations for the entire range of ellipticity ratio examined 
(0,117 < k < 3.7), The apparent advantage in accuracy of the model 
developed by the multilinear regression analysis (Equation 6) com- 
pared with the Hamrock and Dowson equation (Equation 1) (one percent 
versus 30 percent on average) may not justify reconmending Equation 
(6) over the well-known Hamrock and Dowson equation (Equation 1) . 

It is also concluded from this study that the Hamrock and Dowson 
equation can be applied to contacts with the major axis aligned in 
the direction of motion. With respect to the definition of ellip- 
ticity ratio, given by Hamrock and Dowson as the ratio of the 
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semi-major axis to the semi-mino* axis of the contact ellipse, the 
authors would like to suggest that the variable k in Equation 1 be 
redefined as the ratio of the axis of the contact ellipse perpen- 
dicular to the direction of motion to the axis parallel with the 
direction of motion since the former definition precludes a value 
of k less than one. This would permit the application of the 
Hamrock and Dowson equation for contacts with the major axis of 
the contact ellipse aligned in the direction of motion, an orienta- 
tion which has important implications but which is outside the stated 
range of applicability for their equation. 
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IV. SIMULTANEOUS TEMPERAFURE MAPPING AND TRACTION 
MEASUREMENTS IN EHD COOTACTS 

The research presented in this section was in part conducted by 
Mr. E. H. Kool as his Master's degree thesis research [22], It 
was also presented at the 1979 Leeds -Lyon Conference and will 
be published in the conference proceedings in substantially the 
present form. 


Introduction 

One perspective of the contact temperature problem is to divide 
the tenperature rise in the contact into two parts : 1) the bulk 

temperature rise, and 2) tJie local flash temperature rise. The 
bulk tenperature is "a temperature representative of a fairly 
uniform level of those parts of the temperature fields in the 
rubbing bodies that do not lie too close to the conjunction zone" 
[island the flash temperature is "the maximum temperature (rise 
above the bulk temperature) occurring in the conjunction zone" 

[18]. Past research has been primarily directed toward the flash 
temperature component. However, when one gets to specific cases 
both the bulk and flash temperatures are difficult to predict in 
a tribological system. In his 1969 review, Blok [18] presented 
an adaption of the network theory in heat transfer to predict 
bulk temperatures. Since that tune several routine programs have 
become available wliich can be used to predict bulk temperatures 
in complex systems subject to the extent of understanding of 
input parameters such as convection coefficients which are not well 
defined in most tribological systems. 

A major objective of tlie research reported here was to develop a 
technique in which both bulk and flash temperatures could be 
measured in a tribological system as a function of time and space. 
A scanning infrared photometer was employed which covered a region 
including the conjunction region. Several commercial instruments 
were evaluated and are compared. One device was further employed 
in several tribological experiments reported. 


Infrared TemiJerature Measuring Technique 

In the previous infrared temperature measuring work in this laboratory 
[19, 20, 21] a single spot microdetector (Barnes RM2A) was used. 

To automatically scan the region of interest a different IR system 
was required. Several were commercially available and considered. 
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A brief discussion of tlie relative diaracteristics o£ several 
instruments is presented in ^jpendix C. The discussion belos* 
pertains to the AGA 750 whidi was used primarily because it was 
available to us at no cost from the Sdiool of Ardiitecture at 
Georgia Tedi and met most of our requirements. 


Radiation Contributions 

Figure 25 shoi^s the various sources contributing to the radiation 
received by the detector from tlie scanned area of tlie ball in the 
El ID shnulator. 

The four contributions N£, and are tlie nonattenuated 

values of radiations emitted by the ball, fluid, sapphire and 
the ambient respectively. Eadi of tlie above contributions is 
associated witli aai attenuation factor which includes absorjJtion 
losses in varioas media and Fresnel reflection losses at inter- 
faces between any two media. 

The ball surface, being opaque, only emits and absorbs. The fluid 
and sapphire are partially transparent and partially absorbing 
media. The aiiibient radiation refers to tlie background radiation 
from surrounding objects in the room reflected by the ball. The 
radiation collected by the detector- lens is the sum of all the 
above contributions and is given by 

N - ti^N^ + npNp * + t,„N„ . 

One infrared filter was used to partially eliminate the contribution 
from tlie fluid, and is chosen based on the spectral emission charac- 
teristic of tlie lubricant used (Fluid Nl, See y^pendix I) for lubri- 
cant definition). Tlie spectral cliaracteristics of the lubricant 
at 20C and tlie IR filter are shown in Figure 26. The emission curve 
for tlie fluid corresponds to a film thiclmess of 25 ym sandiviched 
between two 16 mm thick sapphire disks. This unit shows that the 
lubricant emission is in the range of 3.1 - 3.7 ym with peak emis- 
sion at 3.4 ym. By using the wide band, high pass filter, all IR 
emission below 3.75 ym, thus all fluid omission is eliminated. 

The dotted lines in Figure 26 show the monochranatic black body 
radiation as a function of wavelength at two different temperatures , 
plotted in arbitrary units. These curves show that the spectral 
region of interest in tliis study is to the left of the peak in 
black body radiation, even at temperatures of 177C. The ball radi- 
ation, which can be considered grey body radiation (emissivity is 
constant over all wavelengths) is therefore more predominant in the 
region of 3.75 - 5.5 ym compared with tlie region below 3.75 ym. 
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Radiation Analysis 

The ball radiation is given by the expression Ngg (Tg) . 

The quantity is the emissivity of the ball surface and has been 

determined to be essentially constant and equal to 0.21. Tjj is the 

ball surface temperature, k Ngg (Tg) stands for the observed black 

body radiation at tenperature T^ after filtering by the IR filter. 

The factor k is a dummy factor to make distinction between filtered 
and unfiltered radiation. Since the variable ^Nbb (T) is dependent 

on the geometry and sensitivity of the detector and also on the IR 
filter characteristics, a calibration was necessary. By using a 
calibrated black body source, the calibration of the modified AGA 
750 unit with three extension rings and the wide band filter was per- 
formed. 

The ball attenuation rij^ is dependent on the optical properties of 

the sapphire and oil. The attenuation factor is determined by 
reflection losses at the sapphire-oil and sapphire air interfaces 
and absorption in the sapphire and lubricant film. Therefore 

Since the IR filter used eliminates all radiation of wavelength lower 
than 3.75 ym, all radiation absorbed in the oil would be eliminated 
at lower tenperatures . This makes Tp(h,Tp) equal to 1. From [11] 

the following data for sapphire transmissivity and fresnel reflec- 
tions are obtained 

Tg = 0.97, = 0.076, P2 “ 0.0064 . 

As total result 

It should be noted that in this situation is equivalent to the 
attenuation of the sapphire radiation (Pg) . 

By using the wide band IR filter, the radiation coming from the 
fluid, Np, is eliminated at low ten?)eratures . Therefore pp * Np * 0. 

The contribution from the sapphire disk is difficult to analyze and, 
as shown before, is very small. Neglecting this radiation would not 
effect the accuracy tolerance of the temperature determination. 
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Aninent radiation, that enters tlie infrared camera lens is mostly 
reflected by the sapphire and ball. This radiation can be divided 
into t\<io parts as illustrated in Fi^re 27. The ’’ambient" (1 and 2) 
can be assumed as a uniformly distributed radiation source witli 
emissivity 1 and radiation level according to that of a black body 
at "room tenperature". Because the IR lens is transparent, and the 
detector is cooled inside, less radiation [19] will be emitted to 
tlie ball and sapphire for reflection. Since the ball is very 
reflective, the image of tlie lens is a dark spot on the ball. 
Experimentally it is found tliat the size of tiiis lens image is 
tliree times the maximum contact diameter. Good alignment will 
provide a situation where tlie EHD contact is in the center of the 
dark spot. 

The amount of ambient radiation from 1 and 2 tliat is observed by 
the detector, can be determined in the following way 

k No >= I' % (To) ' k % (V . 

Because tlie ambient has a uniform temperature distribution, it acts 
like a black body at room temperature; the emissivity is therefore 
equal to 1. 

The ajnbient attenuation is dependent on the three different paths 
of the radiation as demonstrated in Figure 25. Tliis results in the 
following formulation of 

'’o ' ”l * ■'s (1 ■ '’l’^ "2 * % ’'^F 

1st reflexion 2nd reflexion 3rd reflexion 

Numerical values for these properties were found in [11]. By using 
the IR filter, Tp(2h, Tp) will equal 1 and tlie otlier property values 

inserted in the formulation will give 

Hq = 0.076 + (0.97)^ * (1 - 0.076)^ * 0.0064 
+ 0.79 * (0.89)^ * 1 

The observed ambient radiation outside the dark spot is therefore 
rio k = 0.7076 k Njg (T„) . 

Inside the dark spot the influence of ajiibient source 3 is inportant. 
Because the transmissivity of the infrared camera lens is not known. 
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no analytical determination of radiation reduction could be made. 
Experimentally it was found that tlie loss of ambient radiation at 
room tenperature after filtering valued approximately 0,4 Isotherm 
Units (ISU) . 


The total observed radiation by 
outside dark spot: 
/- 0.187 * k Ng3 (ip + 


tlie detector is as folloivs: 
0.7076 * k Ngg (Tq) riSU] 


kN i in dark spot: 

« 0.187 * k Ngg (T^) + 0.7076 * k Ngg (T^) - 0.4 [ISU] 


Froni the radiation formula outside the dark spot, the influence 
of tlie sapphire radiation can be obtained. Suppose a situation 
where the ball and ambient have equivalent tenperature 


kN = 0.187 * k Ngg + 0.7076 * k Ngg (T^^) 
* 0.8946 * k Ngg (T^) 


Expected is 

kN = 1 X k Ngg (Tq) . 

The difference is the loss of radiation in the sapphire disk. Part 
of this will be emitted as infrared sapphire radiation and part will 
be transformed into a different type of energy and will be dissipated 
in a different way. 


Ball- surface Temperature Determination 

Under operating conditions the scanning device will provide a picture 
of the ball in the HID region and a small part of the sapphire disk 
holder. Because the temperature cliange of this sapphire holder 
during one test is not more than 4®C, its radiation level can be 
determined. The black spots on the sapphire holder can be con- 
sidered to be on the radiation level of a black body at the tem- 
perature during operation and will be referred to as AN(= !^bb^Tsh^^- 
Relative to this level, the level of all other points in the picture 
can be measured with isotherm level settings. Tliis relative level 
will be referred to as AR. The radiation level of any point on the 
ball is 
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in the light region: 

kN . i KN * AR - 0.187 * k NggCr^^ji) * 0.7076 * k Njj(T^) 
in the dark spot: 

kN » AkN + AR « 0.187 * k + 0.7076 * k Ngg(Tg^) - 0.4. 

Therefore 


^ ^BB^’^ball^ 


AkN + AR - 0.7076 * k + (0.4) 

DTW"^ 


By using the calibration curve, the ball radiation level will give the 
ball tenperature. 

The IR analysis for the Timken- type friction tester described beloiv is 
less conplicated than that for the EHD device because there is no sap- 
phire plate and the magnification used was much less , hence , eliminating 
the reflected black spot. 


The Tribological Devices tOid Method 

Tivo devices were studied with the IR scanning system discussed above. 

One was a sphere on flat EHD simulator designed specifically f'or tlie 
work and the other was a commercially available LFW-1 Timken- t>T>e 
tester. These two devices are shown schematically in Figures 29 and 30, 


EHD Simulator 

The simulator was a ball driving a sapphire disk. The disk is mounted 
vertically in order to meet the infrared camera requirement of hori- 
zontal radiation detection. In this configuration the ball radiation 
in and around the EHD contact can be exposed to the field of view of 
the camera. The ball is driven through a flexible coupling and the 
normal load is applied by air pressure. When the traction force in 
the EHD contact is larger than the friction force in sapphire support 
bearing, rotational energy from the ball is transmitted to the sapphire 
disk. Both the normal and traction force are measured simultaneously 
by the three-dimensional piezo- electric load cell, mounted immediately 
behind the ball support. Rapid response is possible, because the 
natural frequency of this transducer is 8 kHz. 

The tachometers used for the ball and sapphire speed give 20 pulses/rev 
and 10 pulses /rev respectively. The instrumentation is adequate for 
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tlie determination of the average velocities, but instantaneous 
speed changes cannot be detected. 

After the machine was built and tlie measurement equipment installed, 
three major iiiperfections were found. Tlieir origin and implications 
are discussed. 

1. In order to obtain a small slide-roll ratio, the amount 
of friction in the sapphire support ball bearing was 
reduced by eliminating some balls. Out of tlie 15 balls 
in tlie angular contact bearing, nine were excluded and 
a new separating cage was built for the remaining six 
balls. Teflon was used as cage-material, again to 
reduce friction. Hotvever, the retaining cage left some 
room for tangential motion of the balls. The friction 
betiveen the Teflon and the bearing elements was tlierefore 
not constant. Tliis consequently led to an alternation in 
the traction force and sapphire speed during rolling 
experiments . 

2. Tlie sapphire is mounted in a holder, whicli is mounted 
inside tlie inner race of tlie angular contact bearing. 
Unexpected inaccuracies due to tlie macliining operation 
on the sapphire holder caased tlie sapphire surface not 
to be exactly parallel to the rolling track of the balls 
in the outer race of tlie support bearing resulting in a 
cam action. Because of tliis movement the piston in the 
pneumatic- ^linder also moves. Tlie change in air volume 
in the cylinder together witli changing friction force 
between the piston and the cylinder -wall, were respon- • 
sible for an alternating normal load in tlie contract. 

Fran measurements the maximum piston friction force is 
estimated to be 20 Newtons. Under operation conditions 
the average normal load was varied between 95N and 375N 
to create Hertzian pressures from 1.0 to 2.0 GPa. Hence 
the cam action caused from 5 to 20 percent load variation. 

3. Alignment of tlie ball drive shaft witli respect to the 
motor drive shaft was difficult. In addition, the axis 
of the ball shaft was moving because of the ball movement 
induced by the sapphire wobble. The spring elements in 
the flexible coiqiling created an axial force along the 
ballshaft (x direction). Tlierefore, tlie x direction force 
detected was a canbination of traction force due to spin 
in the EHD contact and the axial spring force. 

Only by using the very sensitive triaxial force transducer were these 
discrepancies between the ideal and real situation detected. However, 
most EHD simulators and traction drives will also have similar imper- 
fections . 
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LFty-X standard Test Device 

This device is a comnercially available Timken>type test device for 
which there is an ASBl Standard Method (D-2714-68) to be used for 
calibration. The results quoted below for it were the result of 
following that calibration procedure, A schematic of the system is 
shown in Figure 30 . The IR scanner viewing direction was coaxial 
with axis of rotation of the ring. The tenperatures reported from 
IR method were of tlie side of the test block as close to the con- 
tact as could be discriminated and the lock-nut holding the ring 
on the shaft. 


Experimental Results 

EHD Configuration 

TVo series of measurements were made for two fluids (N1 and Santo- 
trac 50). Because tlie emission spectra in the 2 ym to 5.5 pm wave- 
length range are almost identical, tlie same filter could be used to 
reduce fluid radiation for botli fluids . The fluid properties are 
summarized in Appendix D. 

The temperatures were determined in the following manner (See Figure 31) . 

1. The sapphire holder tenperature is measured with a 
thermocouple at the beginning and end of the experi- 
ment. Usually this difference is less than 3C in 
rolling conditions ar.d no more than 15C in sliding 
conditions. By using a timer, the total experiment 
duration and me time at which data is taken are 
known. The temperature of the sapphire holder at 
the time of observation can be estimated. 

2. Relative to this radiation level, the radiation 
levels in the center of the contact and in the oil 
wake can be determined. At these two places the film 
thickness is smallest, therefore possible fluid 
influences are minimal. 

The radiation level can be set on the scale to the left of the picture, 
causing all points of this radiation level to be bright white. In 
Figure 31a the radiation level is set for the sapphire holder and gives 
a value of 0.05 on the scale. In Figure 31b the radiation level is set 
for the oil wake and we observe this level to be 0.31. "nie difference 
in radiation levels between the sapphire holder and the oil wake is 
equal to the difference on the scale , multiplied by the range ; thus : 
(0.31-0.05) * 20 “ 5.2 isotherm units. In the same way the difference 
between the EHD center and sapphire holder is determined. By using 
the procedure described in the previous section, the absolute tempera- 
ture and the tanperature difference on the ball surface, between tlie 
EHD center and wake , can be obtained. 
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Measurements in RolUnfi Conditions 

Traction and temperature measurements were simultaneoasly taken 
in the low slide-roll ratio range. The operating conditions varied 
in rolling speed from 0.5 m/sec to 1 m/sec; in load from 1,0 GPa to 
1.9 GPa maximum Hertzian pressure. The slide-roll ratio was always 
less than 0.07, 

Becaase the teiif»erature rise was expected to be small (less than 
IOC) , some of tlie measurements were performed witiiout the IR filter. 

In tills way, the total amoinit of radiation received by tlie detector 
would increase by approximately 100 percent and the difference in 
radiation levels could be detemined more readily, 

The maximan error observed in the unfiltered case was 2.25 isotherm 
units (ISU). Assuming a ball tenperature of 40C, the temperature 
error range is then four degrees C. For the highest load and maxi- 
mum speed the radiation difference AV was equal to 1.6 (ISU). Since 
this value is smaller than the error band, no discrete values of 
temperature rise can be given. Nevertheless, it can be stated from 
tlie obsenrations that witii botli fluids the temperature rises under 
the previously mentioned BHD conditions are less than six degrees C. 

From the force measurements tlie following results were obtained. In 
situations closest to pure rolling, no brake used, the average trac- 
tion force in all runs was almost constant and equal to six Newton. 
IVlien the brake was used, slide-roll ratio and traction force increased. 
Tlie maximum values for traction coefficient were witli sane braking 
0.08 for N1 and 0.10 for Santotrac 50. 


Measurements in Sliding Conditions 

Tlie operating conditions in tliis series of measurements varied in 
ball surface speed from 0.5 m/sec to 1 m/sec and in maximum Hertzian 
pressure from 1.2 GPa to 2.0 GPa. Tlie infrared filter was always used. 
Tlie experimental results of temperature difference are given in Figures 
32 and 33. 

The system thermal transient effects required about 300 sec; the data 
reported were taken in steady state thermal conditions. Tlie equilib- 
riiP’* traction coefficients were about 0.1 and 0.05 for tlie Santotrac 
50 and N1 respectively. 

Because tlie radiation levels of EHD center and the reference (sap- 
phire holder) must appear on the same picture, a larger IR range 
setting (Figure 31) must be used with increased temiierature difference. 
The temperature eiror band due to observation was, in tliis case, also 
4C. During the radiation observation, a horseshoe shaped area of 
elevated radiation was found around the EHD contact. Tlie maximum 
levels appeared at the ends of the horseshoe. Figure 34 gives an 
example of tliis observation. 
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An explanation for tliis appearance was found with help of Dr. Lauer 
[23], The onission spectra of fluids are the same as the absorption 
spectra and are dependent on temperature and pressure. By increasing 
the tein'jerature , the emission band around 3.4 pm tends to shift to a 
longer wavelength and also broadens . An approximation for this 
broadening effect is given in [24] . The halfband width is expressed 
by 


^"^1/2 “ + A exp (-Uor/^T) 

residual width 
o 

Uqj^: mean reorientation potential barrier 

At 70C fluid temjoerature , it can be expected that some fluid radi- 
ation will pass the wide band filter. 

In sliding conditions lubricant temperature iust outside the contact 
area will be equivalent to the ball temperature in that region, which 
vary between 60C and 150C. Because the fluid radiation is a function 
of temperature and film thickness , the effects of fluid band broaden- 
ing in the EUD center and wake are small ; but outside the EHD region 
where the film is thick, the fluid radiation becomes observable. The 
accumulation of this hot fluid explains the two maxima in the horse- 
shoe. 


LFW-1 Configuration 

The results presented from this configuration consist of IR surface 
temperatures of the friction block near the contact zone and the lock 
nut holding the ring in place. Tenperatures of the block center and 
oil bath measured by thermocouples are also shown. These and friction 
as a function of time during a standard ASTM calibration test are 
shown in Figure 35. The calibration test is for steel-on- steel at a 
sliding speed of 7.9 m/min (26 :Q)m) , initial Hertz pressure of 372 
MPa (54 kpsi) and a bath tenperature of 43C (IlOF) . The specified 
oil was a "white mineral oil" unformulated. The EHD film sipport 
is negligible and wear extensive under these conditions. The band 
on the block contact temperature during the early part of the test 
is the result of changing emissivity during that time because of 
buildip of oil and wear debris which increases the local emissivity 
from 0.7 of the initial dry block surface to 0.99 of the dirty surface. 

The results show that the conjunction region is always at least 5C 
above the block center as read by the thermocouple. The tenperatures 
are also influenced by both the friction process and the bath 
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teiiperature whicli on tliis particular device cycles through a 13C 
amplitude. 


Discussion of Results 

The small temperature rise in rolling conditions does not, luider the 
most scvei'c circumstances studied, exceed seven degrees C in the 
work reported in [21]. This is in agreement with the result foimd 
in this work. 

However, in sliding conditions the comparison gives a discrepancy. 

The temi^erature- difference values from [19-21] are much higlier tlian 
tliose reported in Figure 32 where tlie same lubricant and operating 
conditions were employed. Two arguments arc used to exiilain tliis 
discrepancy. 

First, it should be noted tliat different detectors were used. The 
diameter of the spatial resolution of the Barnes BM2A (used in 
[19-21]) is 38 pm £uid approxiinatoly 3.3 times siiuiller tlum the modified 
AGA 750 (spacial resolution: 125 pm). The area over whicli the radi- 
ation is integrated to obtain a single -spot output, is therefore 11 
times smaller for tlie RM2A, and tlie spacial resolution of tlie tempera- 
tui’e distribution mudi better. If we consider a t>T>ical highly loaded 
HMD contact in our simulator (pjj = 1.5 GPa) , tlie diameter of the con- 
tact area is approximately 500 pm. Tlie AGA 750 minimum spotsizc is 
about six jiercent of tlie contact while tlie I^12A minimum spotsi zc is 
about 0.5 percent. It was found in the previous work that the 
thermal gradients in the Hertz region can be very large. The 
spotsize area where the maximum teiiperature is measured, is 
therefore bigger with the AGA 750. The integration over this 
area will level the peak radiation nwre. 


Secondly, a closer look at the surroimding conditions shows a differ- 
ence in base temiierature . The temperature-controlled fluid bath in 
[19-21] provided a constant ball and fluid temperature, except in the 
contact. Usually this temiierature was 40C. In the measurements 
reported here, the base temperature was the ball tenperature reached 
after the thermal transient. This tenperature is dependent on the 
operating variables. Tlie range of tliis base temperature varied 
from 60C to 140C. At these high temperatures the traction coeffi- 
cient is lower and less heat generated in shearing the film. Tliere- 
fore a lower temjierature difference is exiiected. 

Based on both calculations and previous measurements the film thickness 
to surface roughness ratio was from 4 to 40 in the EHD measurenrents 
reported here. Tlie difference in traction coefficient between N1 
and Santotrac 50 is clearly reflected in the temperature rise. Tlie 
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higher traction of Santotrac 50 results in a higher contact tem- 
perature rise as would be expected. 

In the LM-1 result we see the influence of tlie system transient 
and tlie importance of the global tliermal characteristics on tlie 
conjunction temperatures. The batli tliermal control forces tlie 
level of the entire systoii to fluctuate. A preliminary analysis 
of tlie system suggests that the dominate mode of heat rejection 
is conduction out the shaft. The entire ring and shaft end are at 
essentially the same tanperature as the block surface very near 
the contact region. 

In general altliougli tlie scanning infrared system permits deter- 
mination of surface temperatures both in tlie contact and the sur- 
rounding region, the metliod of signal conditioning and display 
makes it difficult to follow the rapid transients and steep gradi- 
ents anticipated in tribological systems. Tlie scanned IR si^al 
inherently contains the information for better time and spatial 
resolution but must be handled differently to take full advantage 
of it. We are currently developing an analog recording capabilitv’' 
for post analysis of tlie signal whicli will permit a more complete 
display of tlie infonnation. 


Conclusions 

An infrared scanning system has been adapted to two tribological 
systems. In many respects the scanner expands tlie versatility of 
the IR teiiperature mapping but tlie added complication and cost are 
not readily justified. The results are consistent witli previous 
more detailed IR temperature mapping whidi is more tediously 
acquired. Further development of the signal handling approacli 
to tlie scanning metliod is expected to make it mudi nwre valuable 
to tribological studies. 
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A FILM THICKNESS ANALYSIS FOR LINE CONTACTS UNDER 
IW ROLLING CONDITIONS WITH A NON-NEOTONIAN 
RHEOLOGICAL MODEL 

Introduction 

Tliis .'inalysis was conducted by Mr. Burak Gecim as an extension of 
his Masters' thesis researcli and will be presented at the 1980 
ASNDi/ASLE International Lubrication Conference and published in 
substantially this fonii in the Transactions ASME Journal of Lubri- 
cation Technology. 

High viscosity fluids, especially when subjected to high rolling 
speed conditions may shw a discrepancy in film thickness, from 
the predictions based on tlie Newtonian theor>'. Dyson and Wilson 
[21] proposed a power- low fonu, nonlinear constitutive equation to 
represent the shear rheological behavior of silicone fluids ajid 
compared the theoretical predictions witli their experimental film 
thickness data on these fluids. But, it should be realized that 
lubrication with these high viscosity lubricants, especially under 
low loading (as in their exjieriment) is hardly expected to be in 
the EHD regime, ajid increasing tlie rolling speed results in further 
deviation from the EHD regime on generalized lubrication regime 
configurations , Figure 36. This anfciguity associated with lA'^son 
and Wilson's analysis indicated the need of a generalized analysis 
of fihii til ickness" behavior with fluids which behave predominately 
non-Newtonian over all regimes of operation. 

Tlie present analysis utilized the nonlinear constitutive equation 
proposed by Winer and Bair [25] , which was modified and applied 
in an EID film thickness analysis by Gecim and Winer [28] for low 
viscosity lubricants. In [28] it was shown that low limiting shear 
stress parajneters, which are material properties, and/or high 
sliding speed? are tlie major causes of decrease in film thickness 
from Newtonian theories. In tlie present analysis of high viscosity 
fluids with lower limiting shear stresses, derivations are confined 
to tlie pure rolling case and tlie objective is to see the effect of 
the low limiting shear stresses on film thickness. The analysis 
is carried out for a line contact geometry under pure rolling, in 
fully flooded full film regime, with the assumption of isothermal 
conditions . 

Fluids of high viscosity may have a low limiting shear stress 
because as the material pressure is increased (or temperature 
decreased) at a given rate the diaracteristic (or relaxation) time 
of the high viscosity material will become equal to tlie character- 
istic time of tlie process at a lower pressure (or higher tenperatnre 
than would be the case for a low viscosity fluid) . Tliis will result 
in the liquid-solid transition occurring at a lower density than for 
a comparable low viscosity material. Tlie lower "frozen- in" density 
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will result in a lower limiting shear stress. Although this 
characteristic is generally expected and experimentally observed 
for some high viscosity silicones [ 2 ] , additional measurements 
of limiting shear stress parameters for high viscosity fluids 
should be made. 


h 


Nomenclature 


Equivalent modulus of Elasticity 


= 2 


A - vj 1 - v'\ 


Dimensionless parameter for Figure 36. 




gj Dimensionless parameter for Figure 36 , 


^ (w/L)^ \ 


1/2 


Film thickness 


m 

P 

% 

p' 


Nominal film thickness 


Dimensionless film thickness parameter 





p. 





slope of the limiting shear stress -pressure relation = mp 
Pressure generated in the fluid 


Maximum Hertzian pressure 


= 0.42 


/w/L • E^\ 

r^/ 


1/2 


Pressure gradient in x^ direction 
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r -1 ^ * ^2! 

R Equivalent Radius R - 

Uj^ Lubricant velocity in Xj^ dii'ection 

{i Surface and rolling speed 

w/L Load per unit length of contact 

and Xj Parallel and perpendicular to the plane of film coordinate 
axis respectively 

T Shear stress in the fluid 

Tg Shear stress on the surfaces 

Tj^ Limiting shear stress of fluid 

a Pressure coefficient of viscosity 

y Viscosity 

Zero pressure value of viscosity 
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Shear Rheological Model and Derivations 
of the Governing Equations 
The non-Newtonian constitutive equation 


dU| 

33^ 


y 


tanh 


'(k) 


is written in terms of its Taylor series expansion as 

3 . . 5 


or 


dUj 

cBc, 


y 




+ . . * 


dUi _ Tl ^ 1 h\ 

dxj ” li (2n-iy yhj 


2n-l 


and then coupled with the integrated momentum equation 


( 1 ) 


( 2 ) 


( 2 ’) 


X = X 3 for -h/2 ^ X 3 < h/2 (3) 

where 

T = t(Xj^,X3) = 0 at X3 = 0 

as shown in Figure 37. Then the velocity distribution predicted 
by the model, is obtained by integrating Equation (2’) subject to 
Equation (3) and is 
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(5) 


2n-l r , > 2n 

,2n /h 


] 


(4) 


vdiich is integrated in the continuity equation 


+h/2 


-h/2 


u, dx- ® u h 
1 o m 


(5) 


dP 


where when = 0. Then the governing non-Newtonian equation, 

corresponding to Reynold's equation of the Newtonian model is found as, 


1 1 /PVhy*^'^ ^ ^ 


E l 1 /P' 

m:2n?T 
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( 6 ) 


or 


tanh 


•1 /Mi\ 


1 1 

7 /FIT 




tanh 




_ 2}ju 
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( 6 ’) 


dP 


tP’h, 


where P' = and ig = |■“ 2“1 is the shear stress on either surface. 

The surface shear stresses are equal in magnitude but opposite in 
sign under pure rolling with the configuration shown in Figure 37. 

T 

If I — 1 = 1 — I is about 0.6 or less, then the LHS of (6') can 

max ”l o 

be approximated by ■=• | — | which yields the classical Reynolds* equa- 

O Tt 


tion. For 0*6 < |— | 5 1, the LHS of (6*) can be approximated by, 

L 
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1.616' 

0.946 

'^S 



''l 



within 151% Error. 


With this approximation the governing equation will take the form 



4.88(pu)^ 


(h-hj' 

m 


h 




(7) 


where 


C 


1 


0.62 . 


It should be realized that in entering the converging wedge of a 
contact, the shear buildup in the fluid passes through stages where 

it is still Newtonian with j— j< 0.6, but for tlie purpose of sira- 

plifying the formulation it is assuned that the fluid is predomi- 

Tg . 

nantly non -Newtonian in the contact entrance, with | — | - 0.6. 

Obviously high viscosities and high rolling speeds siq^port this 
assumption. Newtonian-non-Newtonian transition in terms of the shear 
stress to limiting shear stress ratio is a, more complex phenomena, 
Since at each x^^ location 


T(Xg) 


< 


Ts(Xi) 


with rCXg) = 0 at Xg = 0. Hence, for tlie sake of simplicity, equa- 
tion (7) is considered to be the governing equation. 


In the following section the dimensionless film thickness equation 
for those four different regimes are presented with a given set of 
physical operating conditions. The appropriate regime, hence, which 
film thiclmess equation is valid, is determined from the criteria 
given in Hooke [26], with the following exception. In reference [26] 
it is noted that for values of h’ greater than 11.31 (See Figure 36, 

where h' = ) , although there is a localized flattening of 

the film thickness in advance of the pressure spike, the bulk of 
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the contact retains its cylindrical shape and is regarded as being 
essentially rigid. In the present study that region with h’ > 11.31 
is considered to be in EHD region resulting in a larger EHD area 
and a smaller rigid surface-variable viscosity area. 

In fact, the present R-V equation being a multiple of the R-I 
equation, even in the neiv borders defined as above, is accurate 
away from the transition lines between the regions. Hie R-V region, 
being relatively more coiqplicated to formulate but practically less 
inportant than tlie other regimes, is believed to be fairly approxi- 
mated. This point will further be explored in the discussion. 


Film Tliickness Equations 

Solyion of the Govepiing Equation 
in EriP Regime_with Elastically Defomed 
Surfaces and Pressure Itependent 
Viscosity (E-V) 

In this lubrication regime viscosity changes witli pressure as 


and the limiting shear stress is expressed as 

= mP , 

where the zero pressure value is dropped since it is relatively less 
important than the slope m in expressing the pressure dependency of 
the limiting shear stress, especially in the high pressure region 
of practical interest in lubricated contacts. 

Solution of the governing equation, Equation (7) with these expres- 
sions is outlined in Appendix E and the resulting dimensionless 
film thickness is 
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where 


I(P) « 



g-0.62aP 


and 


I(P) 

a 


-1.05 


0.5 


The numerical value of 0.5 for the bracketed term results from the 
use of an inlet pressure (P^) of 0.68 GPa and a range of repre- 
sentative values of a. This numerical value is quite insensitive 
to the values of both a and P^^ in the ranges of tliose variables 

usually associated with EHD. 


Solution of the Governing Equation in the Elastically 
tteformed Surface "Isoviscoas Regime (E-I) 

In this isoviscous case, y = Uq and the dimensionless film thick- 
ness equation is 





( 9 ) 


since 


I(P) 


p 0.62 

"i 

“or 


for this case. It should be remembered that this regime is the 
lubrication of soft materials with relatively low elastic moduli, 
such that the pressures generated are low enough not to cause an 
increase in viscosity but do cause elastic deformation. A limiting 

1 oP 

case of maximum P. might be in the neighborhood of — so that (e ) 

X Otr 

does not contribute an order of magnitude change to the viscosity. 


Solution of the Governing Equation in the 
Rigid Body Isoviscous Regime. 

In tliis regime pressures generated in the contact are considerably 


I 
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lower, so that neither elastic deformation nor significant change 
in viscosity will occur. As outlined in Appendix ( F ) the dimen- 
sionless film thickness equation for this regime is 



1.738 


„0,469 

m 


0.765 



( 10 ) 


s 

I 

! 

• Solution of the Governing Equation in the 

I Rigid Body - Variable Viscosity I^egime. ' (R-V) ' 

' As in the Newtonian case and outlined in Appendix ( G ) , the film 

i thickness in this regime is expected to be a multiple of the film 

i thickness in rigid body-isoviscous case. The inlet pressure values 

I of interest will be larger than - which will cause significant 

I ® ! 

f changes in viscosity but are bounded to a range which will not cause I 

significant elastic defoimation. This range of inlet pressure can ^ 

1 

( be of the order of 1.0 to 5.0 times -. Therefore for a range of > 

- 1 - 4 

pressure viscosity coefficient of from 14 to 45 GPa (1 to 3 x 10 
“ 1 

psi ) the constant multiplier ranges from 2 to 3. Hence, for all 
practical purposes, this constant is taken to be 2.5 and the dimen- 
sionless film thictoess equation for this case is 



Rigid body-isoviscous 


(11) j 


The film thickness equations presented above are plotted in Figures 
38 through 43 within the ranges as shown in Tables 13 through 16. 
Viscosities with three different orders of magnitude and two loads 
with one order of magnitude difference are used. For each case 
rolling speed is varied within two orders of magnitude , and the 
limiting shear stress parameter m is assigned three different values. 

The corresponding Newtonian film thickness equations for each regime 
are listed for conf)arison [29] 
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Note, for these figures the contact materials were taken to be steel 
on steel with = 220 GPa (= 3.3 x 10^ psi). Therefore, elastic- 
isoviscous regime is not involved in the exanple plots. 


Generalized Regime Charts 

By defining gj^ and g^ as defined in Reference [26], and by letting 



the four different film thickness equations can be represented as 


= 1.738 

R-I 

(13a) 

• 2.5 » h^_j 

R-V 

(13b) 

= 1.24 n,"-'* 

E-V 

(13c) 

= 1.5 C g»-S= 

E-I 

(13d) 
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/E /E \ 

where C « I pf 1 and I I is assimed to be in order of one 
to ten for these E-I conditions. 

Note, however, that the Figures 44 and 45 are only approximate since 
in deriving Equation (13c) above the power of E^ and a are not 

satisfied by the definition of and g^; likewise, in Equation (13d) 

powers of (p^u) and w/L are different then the ones in the original 

^nation. Hence, the purpose of including these two charts of regimes 
is just to get a rough idea of which lubrication regime is to be 
expected under a given set of operating conditions. 

Finally it should be noted that the regime charts for the present 
non-Newtonian model have the same characteristics as the charts 
based on Newtonian models [26], except that the extension of E-V 
region, and therefore, R-V region being confined to a smaller area. 
This point will further be ejqilored in the discussion. 


Discussion 

The shear constitutive equation proposed by Winer and Bair [25], 
has an elastic term and a nonlinear viscoplastic term to describe 
the shear rheological response of the fluids in EHD contacts, that 
is 


Y = 


1 dt 


at ^ 



High viscosity fluids, with relatively low values are expected 
to have some elastic characteristics of behavior, therefore, it was 
first attempted to solve the full model in a conventional film thick- 
ness analysis. 


dr dr 

The term ^ can be written as where ^ and since 




Xj + F(Xj^) from the momentum equation 
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and therefoie 


dr ~ u dt _ u 

at ^ a5T " ^ 
00 00 00 


d^P 

— x-X, 

dx^ ^ 


dF(x^)" 

dxj^ 


Inclusion of the second derivative of pressure in the model results 
in extremely high pressures (and pressure gradients) in the inlet 
zone and therefore causes numerical problems. Although the approxi- 
dx, 

mation = u can give some idea for conparing the relative effects 

of elastic and viscoplastic terms of the constitutive equation, the 
velocity profile is not constant across the film and this assurrption 
may be causing the numerical problems when the whole constitutive 
equation is twice integrated across the film. 

The apparent viscosity versus shear rate curves in Reference [27] 
imply a visco-plastic behavior for high viscosity fluids (although 
the viscous region is confined to very low shear rates) , and resem- 
bles the limiting shear stress -plastic deformation characteristic 
of visco-plastic behavior as presented in [25], by Winer and Bair. 
Therefore it is assumed that the visco-plastic portion of the shear 
constitutive equation. Equation (1), can be used in the film thick- 
ness analysis with high viscosity lubricants. 


The governing equation. Equation (6*), is approximated by a simpler 
function within +5 percent error in order to make the formulation 

T T 

easier. This ^proximation is valid for | — | ^ 0.6. For | — | <0.6 


behavior is purely viscous and the governing equation is the well- 
known Reynold's equation. This fact is also seen from the shear 


less than 


constitutive equation because for the argument 

s approximately e^al to |;^| and 


approximately 0.6, tanh’ 


T 

Tr 


the constitutive equation is simply Newton's law of viscosity. The 
non-Newtonian characteristic of the proposed shear constitutive model 


will be dominant for [™( 5: 0.6. It should also be noticed that both 


the velocity distribution. Equation (4), and the governing equation. 
Equation (6) , can be reduced to their Newtonian counterparts by 
t^ing the limit approaching infinity. As explained above for 

tliis characteristic, |;^| < 0.6. is sufficient. 


After deriving the governing equation by following the conventional 
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procedure of continuun mechanics, this equation is solved for the 
filjn thickness for several cases. These are, as outlined in the 
introduction, tlie four different regimes of lubrication. The 
inportance of discriaiinating between these regimes is seen by 
noting that the difference between Martin's solution and Dowson 
and Higginson's solution will be more tlian 50 percent, for the 
data of Figure 42. This difference, however, will decrease as the 
regime shifts from R-I to R-V and subsequently approaching E-V 
region. 

Early studies of R-V regime indicated that the film thickness in 
this regime can be expressed as a multiple of tlie film thickness 
of R-I regime, and the multiplication factor is a function of the 
pressure coefficient of viscosity [30]. However, in Reference [30] 
it is pointed out that the limit of this factor is approximately 
2.3. Beyond this limit, the solution for this region is referenced 
to some extrapolations of Blok's results [31]. However, in Green- 
wood's paper [31], it is shown tliat the extrapolation will not 
deviate mucli from the Grubin's E-V film thickness equation. Also 
in Reference [26] it is noted that for the cases which result in 
factors larger than 2.3 there is a localized flattening of the 
film in advance of the pressure spike. These all inply the possi- 
bility of extending tlie E-V region, and confining the R-V region 
to a smaller extent, as presented in this analysis. Hence, for 
the non-Newtonian R-V region, whidi is believed to be of relatively 
less practical importance and more difficult to formulate, tlie film 
thickness is assumed to be I S times the film tliickness of R-I 
region. 

The film thickness equation for E-I region is not used in the 
numerical exanples presented because the contact materials are 
considered to be made of steel ana therefore, the pressure generated 
in the film cannot be high enough to deform steel and yet low 
enough not to cause significant cliange in viscosity. In drawing 
Figures 44 and 45, the value for E^/P^ is assumed to be in the order 

of 1 to 10 because of the low modulus materials involved. Tlie 
border lines of the regimes are only weakly dependent on (Ej./Pj^) 

of E-I case. As explained in the introduction the limiting shear 
stress parameters tend to decrease with increasing viscosity. 
Although this characteristic is generally accepted the present 
analysis suffers from the lack of experimental data for the slope 
m values in tlie high viscosity ranges used. In a previous study 

-1 -2 -"i 

[28], for the viscosities in the order of 10 to 10 Pas [10 

or 10 Ibs/in ], the experimentally measured slopes of 0.1 or 0.05 
have been used and it has been observed that the zero pressure 
value of T, is relatively less important tlian the slope m. In the 

Xj /I 9 

present study with viscosities ranging from 1 Pas [10*^ Ibs/in^] 
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to 10^ Pas [lO"^ Ibs/in^] , slones ranging from 0.01 to 0.001 are 
thought to be reasonable. In [32] capillary viscosity measurements 

2 

indicate that dimetliylsiloxane (10 Pas at 27C) has a limiting 
shear stress of 4 MPa at 550 MPa confimiLng that tlie slope is in 

the order of 10 . Tliis measurement was at low temperature and 

increasing tlie teuperature is expected to result in furtlier decre- 
ment of the limiting shear stress. Finally, it should once again 
be rememl^ered tliat t!ie order of the slope m is very inportant in 
predicting tlie film tliickness since tlie dimensionless film tliick- 
ness is directly related to viscosity and m, tlie effect of increase 
of tlie former can only be compensated, other parameters remaining 
the same, by a decrease in the latter. Hence, substantially lower 
film thickness measurements with high viscosities imply the neces- 
sity of assigning low limiting shear stress parameters in such 
analytical studies. 


Conclusions 

In a recent paper Gecim and Winer analyzed tlie effect of lubricant 
limiting shear stress on film thickness, based on a non-Newtonian 
rheological model [28]. In that study the viscosity range was 
fairly low so tliat, the deviation (»40 percent) from the Gnibin’s 
conventional tlieory is attributed to high sliding and, to some 
extent, low limiting shear stress paraneters (whicli have been 
measured £uid reported in [2 ]). No drastic deviation from Newtonian 
solution, or a sudden collapse is expected for that case, as con- 
cluded in [28]. However, as the experimental studies indicated, 

EHD contacts failed to form a lubricant film witli high viscosity 
lubricants and tlie results of the present analysis imply this 
analytically. Since tlie theoretical film thickness equation of 
the present analysis is directly related to the viscosity, it is 
concluded that, for this pure rolling case, these deviations from 
the Newtonian theory are due to their low limiting shear stresses. 

For eadi viscosity case, with low loads where the regime falls into 
R-I region, greater deviation is predicted than with high loads 
where the regime will be in E-V region. It should be recalled that 
the R-I region implies low contact pressures, hence, low limiting 
shear stresses. Increase of rolling speed for each curve results 
in an increase in tlie deviation of non-Newtonian lines from tlie 
corresponding Newtonian line. 

Thermal effects, which might be considerable with high viscosities, 
indicate the need of inclusion of the energy equation in the non- 
linear formulations of the present analysis. The authors do not 
address that question in the work. 
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Although the present isothermal analysis is not able to match 
the experimental data of Reference [27] with the m values used, 
it should be noted as it is reported in Reference [27] that the 
test fluids with these viscosities were too viscous to circu- 
late into the contact and, some fracture was observed at high 
rolling speeds. This ten^ to confirm the inability of these 
high viscosity fluids to form a satisfactory lubricant film, due 
to their low limiting shear stresses which might be readied even 
far out in the inlet zone. 

Finally, it is concluded that, with steel on steel line contacts 
the lubrication with high viscosity fluids falls into R-I regime 

under the loads of the order of 10^ N/m [10^ Ib/in], and into E-V 

region if the loading is in the order of 10 N/m [10 Ib/in] or 
higher for the rolling speeds of practical interest. The inter- 
mediate R-V regime is confined to a smaller region and thought 
as being less important. 

The present study supports t!ie observation of film thicknesses 
which are smaller than predicted by traditional lubrication analysis 
for high viscosity fluids and implies the need for 

i) including the energy equation in the nonlinear 
foimulations , and 

ii) measurements of limiting shear stress parameters 
for high viscosity fluids. 
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Appendix A 

Experimental Fluids 

ejection 11) 


SYMBOL: 

SOURCE: 

TYPE: 

PROPERTIES: 


R620-15 (See also Appendix D) 

Sun Oil Company 
Naphthenic Base Oil 

2 

Kinematic Viscosity at 37. 8C m /s 
24 X 10"^ 

2 

Kinematic Viscosity at 98. 9C m /s 
3.7 X 10'^ 

Density at 20C kg/m^ 

0.916 X 10^ 

Average Molecular Weight - 305 


SYMBOL: 

SOURCE: 

TYPE: 

PROPERTIES: 


R620-16 (See also Appendix B) 

Sun Oil Company 
Naphthenic Base Oil 

2 

Kinematic Viscosity at 37. 8C m /s 
114 X 10"^ 

- 2 

Kinematic Viscosity at 98. 9C m /s 

8.1 X 10"^ 

Density at 20C kg/m^ 

0.930 X 10^ 


Average Molecular Weight - 357 
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S'^MBOL: PL4520, PL4521, PL4523 

SCXJRCE: Rohm and Haas Con^any 

TYPE: Polyalkylmethacrylate 

(Polymer additive used in solution in R620-15, 

49% polymer by weight) 

The chemical conposition of each is the same. They 
differ only in molecular weight and are si^plied in 
a carrier oil similar to R620-15 


PROPERTIES: PL4520 PL4521 


PL4523 


POLYMER m: 42.6 


36.11 19.0% 


VISCOSITY 355,000 560,000 2,000,000 

AVERAGE 

MOLECULAR 

WEIGHT 
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Appendix B 

Properties of Steel Rollers and Sapphire Disk 
tSiction III) 


A- 2 TOOL STEEL ROLLERS 

MODULUS OF ELASTICITY E * 207 GPa (30 x 10^ psi) 

POISSON'S RATIO V = 0.3 

HARDNESS OF ROLLER - 59 ROCKWELL C 40 MIN 0 lOOOC (1830 F) 
SURFACE ROUGHNESS = 0.13 vm AA 

SYNTHETIC SAPPHIRE AJI 2 O 3 DISK 
MODULUS OF ELASTICITY E = 365 GPa (53 x 10^ psi) 

POISSON'S RATIO V = 0.25 

HARDNESS Moh 9, Knoop micro indenter 2000 


(=1 63 Rockwell C) 

TENSILE STRENGTH 0.40 GPa (58 x 10^ psi) 

COMPRESSIVE STRENGTH 20.7 GPa (30 x 10*'’ psi) 

REFRACTIVE INDEX n = 1.76 

SURFACE ROUGHNESS = 0.00635 \m AA 

SAPPHIRE DISK - STEEL ROLLER COMBINATION 
EFFECTIVE MODULUS OF ELASTICITY 
E' =286 GPa (41.4 x 10^ psi) 


WHERE 


E' = 2 


-1 


1 - v: 


STEEL 


1 - u: 


SAPPHIRE 


^S'TEEL 


"SAPPHIRE 


COMPOSITE SURFACE ROUGHNESS a = = 0.13 ym 
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Properties of the Lubricant 

— — (S'earon ni) 

NAME: SUNOCX) NAPHIMENIC OIL-2 

CX3DE: R-620-16 

AVERAGE MOLECULAR WEIGHT: 357 

REFRACTIVE INDEX 1.5173 

DENSITY @ 20°C 930,3 Kg/m^ 

0.0336 Ibm/in^ 


VISCOSITY INFORMATION 
IBIPERATURE VISCOSITY 


°C F 

mPas 


REYN (x V 

24.0 75.0 

316 


45.9 

27.2 80.9 

246 


35.7 

40.0 104.0 

93 


13.5 

98.9 210.0 

7.4 


1.1 

PRESSURE VISCOSITY COEFFICIENT 



TEMPERATURE 


a 


°C 

GPa’^ 


Psi'^ (x : 

24 

36.7 


2.51 

27 

35.7 


2.459 

40 

31.9 


2.178 

99 

19.8 


1.365 


TEMPERATURE VISCOSITY COEFFICIENT 

20 < T (°C) <40 6 = 0.047 F'^ = 0.086 C"^ 
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Appendix C 

Conparison of Several Cwnmercial TR Detectors 
[S'ectloii IV3 

The simplest IR microscopic detector is the single spot detector. 
This type, whidi is used in previous work cf. [19,20,21], has no 
automatic scanning option. The integrated radiation coming from 
a spot, 0.001 incli in diameter, is measured. Since the contact 
diameter and region of interest is much larger, an automatic scan- 
ning device is required. 

In order to compare tlie different available devices, some labora- 
tory tests were performed while the instruments were on loan. 

The basic test was to find the temperature distribution on the 
ball-side of a sapphire-steel EHD contact. Necessary data to make 
a comparison is given in Table 12 , Some devices have the option 
of scanning one line in the field of view separately. The renewal 
of this line is usually faster than the renewal of the whole view. 
In case of rapid transient effects, this could be a very useful 
option. 

By analyzing the different features of the detectors, the follow- 
ing criteria could be set to obtain the desirable temperature 
information from a EHD contact. 

a. From the geometrical point of view, a spacial resolution 
of less than 1/1000 indi and a focusing distance of more 
than 1 inch are needed. 

b . The temperature resolution for a gray body (e * 0,2) 
should be better than 1“C (•> Black body resolution 
better than 2®C). 

c. For transient measurement a rapid scanning feature 
would be needed. 

d. For storage reasons a digital output is preferable to 
a picture; combination of both would provide a good 
information. 

These criteria could not be fulfilled by any available detector. 
Comparing the data information, the UTI - 900 unit comes closest 
to the requirements. 
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SYMBOL: 

SOURCE; 

TYPE: 

PROPERTIES: 


SYMBOL; 

SOURCE; 

TYPE: 

PROPERTIES: 


Appendix D 

Lubricants Employed " EHD Simulat or 
(Section IV) 

N1 

Sun Oil Conpany 

Naphthenic base oil R-620-15A 

2 

Viscosity at 37. 8C m /sec 
24.1 * 10‘^ 

2 

Viscosity at 98. 9C m /sec 
3.73 * 10"® 

Density at 20C kg/m^ 915. 7 
Pressure Viscosity Coefficient 
Pressure Viscosity Coefficient 
Pressure Viscosity Coefficient 

Santotrac 50 

Monsanto Company 

Synthetic Cycloaliphatic Hydrocarbon Traction Fluid 

2 

Viscosity at 37. 8C m /sec 
34 * 10"^ 

2 

Viscosity at 98, 9C m /sec 
5.6 * 10"^ 


Density at 37. 8C kg/m^ 889 





Pressure Viscosity Coefficient 

at 

20C 

GPa"^ 

39 

Pressure Viscosity Coefficient 

at 

40C 

GPa“^ 

26 

Pressure Viscosity Coefficient 

at 

70C 

GPa‘^ 

16.7 

Pressure Viscosity Coefficient 

at 

lOOC GPa"^ 

^ 12.; 


at 25C GPa"^ 26 
at 40C GPa'^ 22.5 
at lOOC GPa’^ 15 
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Appendix E 

Derivation of the Dimensionless Film Thickness Equation 
in Elastically Deformed Surfaces , “ 

Variable Viscosity' Regime (i.e. EhD) 

-"TSicHOTiTf 

The governing non-Newtonian pressure gradient equation, Equation 
(7) was 


dP 


“4.88 (y u) 


^ 1 r (b-hj^ , 

' ' ' -rscTT- ‘="r^ ' "-*2 


aP 

Defining y, and Xj^ as functions of pressure, y “ y^e and x^ “ mP 

and assuming that pressure is zero at a large distance from tlie inlet 
point. Equation (A.l) can be integrated 


e_-.^Ld!L=4.88(P Gf ( 

K ° 3.„ 


(h-b^) 

w- 


h 


where bjj^ = h^ the nominal gap of the EHD region. The elastically 
deformed inlet zone film shape 



■/ \2 ^ 

1/2 

- 

/ V 

l/2“|j 

fll 

■ 1 

- iJ-n 

1-1 + 
1 a| 

§-'\ 

Jj 


h ' >’o * k 

can be approximated by the expression 


2 / X 

' ''o * ?ir r 3 )r I ^ ‘ Hr ^ ^ 


as in Reference [27]. 

If we let the LHS of Equation (A. 2) to be I(p^), then Equation (A. 2) 
will take the form 


0 / 2 ^ 3 / 2 \^ 

- c f ( IT ^”3 — ) 

I(PJ = 4.88 Cu u)‘^ a \ V5 — U 

r a^ 2^/2 V2l2cH 


If we change variable of integration on the RHS by letting 


,2 23/2 

^ "T- 


1^/^ = tan^e 


I(Pp-4.S8(P<, kJ 


tan^^6[l + tan^6]tan 9 
[tan^e + tan^-^e 


where 


4 l^^o 3 


The integral c . the RHS is sinply 


fsine^c* 

J (tan^e + l)^'^ J 


1/3 2c - 1/3 
' cos e 


Therefore the governing equation will be 

4/^0 3 \ 

I(P.) = 4.88(y^u) J 


sine ^^■^^/^cose '^/^de 


1 


1 




where tlie integral at tlie RHS, 1(G), is * 0.41 for c = 0.62. 
Therefore 


3 ■ ^173 


pzr 
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recall that c = 0.62, 
Therefore 



o 


Equation (A. 4) is the dimensionless film thickness equation based on the 
non-Newtonian constitutive equation. 
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i^ppendix F 

Derivation of the Dimensionless Film ^tickness 
fequation in Rigid Bocfy Isoviscous Reginte 
XSicH^V)"" — 


The governing equation, Equation ( 7 ) witJi M = will take the form 



2 X 

where ^ * hence ^‘jjj “ ^ ^ ' 

cos y 


If ive lot tanS “ — ^ — , 


I 


TOT = 4.88 u^' 




d0 


where 


tan^ 0 


m 


2 

^m 

mr 


since P = 0 at 0 = - j (i.e. , x = -“) 


and 

P = 0 at 0 = 0 , 1 

m 

[Reynold’s 

the 

value of 0 is found by trial 
in 


0 ... 




“ 


( 

cos^ 0 


•' - 7t/2 



B.C) 

and error until 

4 

cos^ 0 
cos^0 

m 


(A. 5a) 


(A. 5b) 


d0 = 0 
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This value is 


= 0.4525 


then 


P = 


c -c 

•Z 0 7 11.1 0 — y ■ 


0 

_ w, 

c 



COS‘^6 


O.U3 m Po 

3 

cos% 

d0 

■t 


- V2 

L 

- 


^).62 


Since 


m 


W/L = 


then 


d 6 

P TT- , and recall that c = 0.62, 

cos^e 


- 7T/2 


h 1/c 

W/L Po ” , l/c+2 ^ “V/Mq ) 


O 


! 


m 


- tt/2 


J 


- tt/2 


2fi 

cos 0 


4 

cos^ 0 


2T" 

cos 6^ 
m 


1/c 


The value o£ the integral at the RHS is found to be 0.142 with 
1000 grid points of numerical integration between 0 = - it/ 2 and 
0 = 0^ = 0.4525. Then, finally, the dimensionless film thickness 

equation for this case is 



1.738 


0.469 


m 



0.765 


(A. 6) 


64 


Appendix G 

Derivation of Dimensionless Film Thickness 
Equation in kigid Body - Variable Viscosity feecime 
(Section V) — ^ 


In this case, the RHS of the Equation (A. 5b) is the same but the 
LHS pressure integral has tJie form 




dP 


Writing the exponential term in Taylor series expansion and inte- 
grating results in 




P 

o 



dP = P^^ 



which can be approximated by the form, within |5| percent error, 

[AP^* ], which ’will lead to a similar solution as in R-I case, 
where the values of A changes with the pressure range in each a 
case, and with a itself. But, as stated in the text, for all 
practical purposes, the variation of the coefficient in front of 
the rigid body-isoviscous solution, between ” 2.0 to ** 3.0, is 
handled by assuming it to be 2,5. 


Hence, the dimensionless film thickness equation for this case is 




Rigid body-isoviscous 


(A. 7) 
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Table 1. 

Pressure Viscosity Coefficients 


I 

Material 

Tenperature/C 

a^/CPa'^ 

a^/GPa" 

;R620-15 

26 

27.4 

27.7 


99 

15.4 

14.8 


227 

12.0 

8.85 

R620-15 + PL4523 

26 

25.5 

25.7 


99 

17.1 

15.0 

j 

227 

16.8 

10.3 

jR620-15 + PL4521 

26 

24.2 

24.9 

i 

1 

99 

15.0 

15.3 

! 

227 

13.8 

9.8 

R620-16 

i 

26 

35.6 

35.8 

* 

t 

f 

99 

19.8 

19.8 


227 

10.8 

10.6 


where 


a 


J 


p— ♦(» 


dP 


L-'p = 0 


-1 


at constant tenperature 


- 

% “ 3p 


T = . .nst. 


from falling body viscometer data. 
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Table 2. Geoinetxy of the Crowned Rollers 
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Table 3. Dimensionless Speed, Load and Corresponding Measured 
Minimm Film Thickness Results for k = 3.7, R /R = 








Table 5. Dimensionless Speed, Load and Corresponding Measured 
Miniraum Film Thickness Results for Ic = 0.958, R /R = 
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Table 6. Dimensionless Speed, Load and Corresponding Measured 

Minimum Film Thicloiess Results for k = 0.651, R /R = 0.54 


T 
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Table 7. Dimensionless Speed, Load and Corresponding Measured 

Minimum Film Thickness Results for k = 0.360, R ^R^ = 0.21 


r 

I 
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Table 10. Film Tliickness Statistics 


Standard 

Mean Deviation 


iln(H, 

'• nan. exp 

- 10.06095 

0.82934 


^nin. exp 

26.27 pin. 




0.667 pm 



ih . -h. |W 

' min min. exp ‘ 

6.682 pin. 

4.069 

yin. 


0.170 pm 

0.103 

ym 

(h . /h . 

nan. exp' nan.' 

1.302 

0.214 


(h . /h . 

r- min. exp' min.' 

1.0095 

0.1398 



(1) calculated by Hamrock and Dowson Model Equation (1) 

(2) calculated by regression model Equation (6) 
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Table 11. Film Tliickness Statistics as a Function of Ellipticity Ratio (k) 


mn.eiqj. 


k 

Mean 

Standard 

Deviation 

3.70 

1.16 

0.07 

2.40 

1.14 

0.10 

0.958 

1.38 

0.09 

0.651 

1.27 

0.11 

0.360 

1.37 

0.12 

0.305 

1.38 

0.10 

0.117 

1.41 

0.45 

All coiiibined 

1.30 

0.11 


Hamrock and Dowson Model Equation (1) . 


Table 12. Data ior Conparing Conimercial IR Detectors 
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Table 13. Load per unit length and corresponding 
Ma-xiiiiuin Hertzian Ihressure 


Case 1 

87. 6 kN/m (500 Ib/in) 

0.5 GPa (72746 psi) 

Case 2 

876 kN/m (5000 Ib/in) 

1.6 GPa (230,043 psi) 


Table 14. Zero pressure value o£ viscosity in each loading 
case, and corresponding limiting shear stress 
parameters 


Dynamic Viscosity 
at 30OC 

Slope of the Limiting Shear 
Stress -pressure Relation 

0.863 Pas 

(l.;:4 X 10'^ Ibs/in^) 

0.01 

11 .68 Pas 

(1.69 X 10^ Ibs/in^) 

0.01 and 0.005 

91.0 Pas 

(1.32 X 10'^ Ibs/in^) j 

0.005 and 0.001 


Table 15. Physical Input Parameters 


Pressure coefficient of viscosity j 

14.5 GPa"^ (1 X 10"^ psi"^) 

Equivalent Radius 1 

0.0127 m (0.5 in) 

Equivalent modulus of Elasticity j 

. - - ... ... ...... . i 

220 GPa (3.3 X 10^ psi) 

Materials 

steel on steel 

Rolling Speed u 

127 on/sec (50 in/sec) to 
1778 cm/ sec (700 in/sec) 
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Figure 1. High Ten^erature Visccmieter 
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Figure 2. High Pressure Viscometer 










r Stress /MPa 






Shear Stress /MPa 
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Figure 7. Flow Charts for Three Polymer Blends at 26C. (Identical 
Base Oil and Three Different Molecular Weight Polyalkv'l' 
methacrylates.) 



ree 









r 




Figure 11. Pressure-Viscosity Isotlierms for (.) R620-15, 
(o) R620-16, (x) R620-15 + PL4523, (A) R620-15 
+ PL4521 



Viscosity/Pas 


10 



Figure 12, Loiv Shear Stress Viscositv^ by T\vo Tediniqucs for 
(O R620-15, (+) R620-16,' (x) R620-15 + PL4523 
at 26C 
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Figure 14. Crowned Roller and Si^jport Assembly 
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Figure 24 


min 


. Predicted Dintens ionless ^.lini^nun Film 'Iliickness vs 
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Figure 29. Schematic layout o£ the EFD simulator 
and attached equipment 



Figure 30. Schematic of LFW-1 test apparatus 
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Figure 31. l^idiation levels on ball :md sapi^hire holder in sliding 
conditions. Iso-radiation levels are white. 
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Figure 32. Temperature difference between EIID center and wake 
in sliding conditions, for N1 



Figure 33. Temperature difference between EHI) center and wake 
in sliding conditions, for Santotrac 50 


J 
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Figure 34. Infrared picture of sliding contact to illastrate the 
"horseshoe" aroimd the lillP contact. 


Oi^crating conditions: 


'hall ■ "-^5 

o„ =1.75 rjM/m^ 
llz 

Fluid: N1 



Figure 35. hnV-1 system temjieratures and friction: (11 Block- 

contact (IR), (2) Block center (thcrmocouiile) , (31 
Oil hath (T/C1 , (41 Shaft lock nut (IR1 


Coefficient of friction 
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FILM thickness PARAMETER 
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Figure 42. Speed Parameter 


versus Film Thickness Parameter 



W/L = 8Z6 kN/m = 91 Pas = 5 x 10® N/m^ 

1. R-I Newtonian (Martin's Solution) 3. m = 0.005 Model R-I 

2. m * 0.01 Model R-I 4. m = 0.001 Model R-I 



FlUM TWCKIIESS FAKAMETER 
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SFEED PARAMETER 


r 

Figure 43. Speed Parameter U = versus Film Thickness parajneter 

[h4] 


W/L » 876 kN/m 


Wg » 91. Pas = 1.6 X 16® N/m^ 


1. E-V Newtonian (Dowson and Higginson) 

2. m = 0.01 Model E-V 

3. m = 0.005 Model E-V 

4. m = 0.001 ^fodel E-V 
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Figure 45. Regin»e Chart (Non -Newtonian) 
m = 0.001 


